





propaganda purposes or in the popular press. Quotation of titles and abstracts 
in publications of purely technical or scientific nature is not prohibited. 


ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 
Band 23 Dezember 1943 Heft 6 


Inhalt: 








Seite | Seite 
Hauptaufsätze. W. Müller: Über die Berech- | H. Karl: PBiegung gekrümmter, dünnwandiger 
nung der Kräfte an einem in der reibungslosen | 00T ER NG ER d re an a: N 
Flüssigkeit beschleunigt bewegten Tragflügel . . 30 H. Blume: Über die Analyse kurzer Kurvenzüge . 246 


H. Neuber: Die Grundgleichungen der elastischen Kleine Mitteilungen. K. Ludwig: Das Auf- 
Stabilität in allgemeinen Koordinaten und ihre heizen einer Wand durch konstante Wärmestrom 
Integration . . .» ee ea dichten . 





HAUPTAUFSÄTZE 


Über die Berechnung der Kräfte an einem in der reibungs- 
losen Flüssigkeit beschleunigt bewegten Tragflügel. 
Von Wilhelm Müller in München. 


Die Arbeit bringt eine Anwendung einiger allgemeinen, vom Verfasser auf- 
gestellten Beziehungen für die Kräfte an einem in der reibungslosen Flüssig- 
keit bewegten Körper auf die beschleunigte Bewegung eines Joukomwski- 
schen Tragflügels. Es werden Formeln für den Einfluß der verschiedenen 
Beschleunigungsarten, sowie der Drehung und Drehungsschwingung auf die 
Kraftwirkung angegeben. Die Rechnung besteht im wesentlichen in der Aus- 
werlung von komplexen Oberflächen, bzw. Randintegralen nach funktionen- 
theoretischen Methoden. 


1. Das allgemeine Formelsystem. 


In einer größeren Arbeit‘) habe ichh mehrere Formeln aufgestellt für die Kräfte an einem 
in der reibungslosen Flüssigkeit mit beliebiger Beschleunigung bewegten Körper. Bezeichnen 
wir die absolute Feldgeschwindigkeit mit dv, die Translationsgeschwindigkeit mit d,, den 
Vektor der Winkelgeschwindigkeit der augenblicklichen Drehung um eine im Körper C feste 
Drehachse mit &, entsprechend der Umfangsgeschwindigkeit u=®»Xr, ferner mit ® das 
Potential der im Raumgebiet T außerhalb C als wirbelfrei vorausgesetzten Strömung, so ge- 
winnen wir für die resultierende Druckkraft den allgemeinen Ausdruck 


ur R . 
Bo=ozz;| Pndf+gv.xX |vxndf+e \ux(wxndf). 5 ME 


Ü € c 


wobei die Integrale über die Oberfläche des Körpers C zu erstrecken sind und n den Einheits- 
vektor der äußeren Normalen bezeichnet. Die beiden letzten Integrale können auch so um- 


geformt werden, daß die im Innern, bzw. am Rande zu denkenden „erzeugenden“ Wirbel bzw. 
die entsprechenden Zirkulationsgrößen auftreten. Wir erhalten dann z. B. 


Bir i ’ 
Be=o (tax) \ Ondf—-0o0,X \rotvdr— o\uXrotvdr. 

Be. : £ h er 

( To To 
wo Jetzt die letzten Integrale über das Raumgebiet T, im Innern des Körpers zu erstrecken sind. 
1) Wilhelm Müller: Zur Theorie der Kräfte bri der beschleunigten Bewegung eines Körpers in der reibungs 
losen Flüssigkeit, Ing.-Arch. Bd. XIV (1944), S. 32 bis 359. An dieser Stelle ist auch die Literatur ziemlich vollständig 

angegeben worden. 
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Wenn ® und rot gleichgerichtet sind wie im ebenen Fall, so wird 
uxXrotd=@X(tXrotd). 
Bezeichnen wir ferner die resultierende gerichtete Zirkulation mit £/’ und den Örtsvektor 


des sog. „Schwerpunktes der Zirkulation“ mit r,, so können wir schreiben 


\rotvdr=Etl; \rXrovdr=r,Xfl.. 
I, To 
Damit erhalten wir dann z. B. im stationären Fall ohne Drehung den bekannten Ausdruck 
für die Magnus- bzw. Auftriebskraft 
RB. n oO v,xEtf 
In derselben Weise ergibt sich das Moment der Kräfte 


o® 





M=o\rx AL df+o\ıx |, XwXndf)]+oe\rx[uxwxXxudf)] 
© MH 
E oo» > ; ” 
o\rx je" df+o@oxX\rxPndf+teovn,X|\Pndf (4). 
C & 
—o\rX(d,Xrotb)dr— o\lrX(uXrotv)dr 
gi T. 


Im ebenen Fall, wo » auf r und u und u aufr senkrecht steht, verschwindet der zweite und 
letzte Ausdruck, und wir erhalten die vereinfachte Formel 


ar i 
Me=o;; \rx®nds+ov,xX|\®nds—o \rx(w,Xrotp)dr 


D. i ; 
f nds+0o0,X \ DPnds—ov,X \ rXrotpdr 


2 , | a SE 
c 

=o \ ex 3 | 

wobei im letzten Integral wieder 


— 08, X \rxXrotvbdr=— 00, X (u, XEN) 
gesetzt werden kann. 

Wie’ich gezeigt habe, sind die gegebenen Formeln von so allgemeiner Bedeutung, daß 
sie sich unmittelbar auf fast alle Sonderfragen, z. B. aus der Theorie des Propellers, der 
Turbinen, sowie auf die Luftschiff- oder Flugzeugrumpfaerodynamik anwenden lassen. Im 
folgenden wollen wir uns im besonderen mit dem ebenen Fall beschäftigen und zunächst im 
Anschluß an die oben angeführte Arbeit, die Ausdrücke für Pc und M. mit Benutzung von 
komplexen Größen so umformen, daß die vorkommenden Integrale mit funktionentheoretischen 
Methoden unmittelbar ausgewertet werden können. 


2. Die komplexe Darstellung im ebenen Fall. 


Wir gehen am besten aus von dem im Raumgebiet 7 enthaltenen Impuls, den wir mit 
bezeichnen. Nach dem Gaußschen Satz ergibt sich, wenn wir die innere Normale «; nennen, 


S=o\dbdt=o \r $dt=o | Pnds=—o\Dnds...... . 
Mit 
u=cos(nz)e,+sin(nz)e,; dscos(nx)=dy, dssininx)=—dx 


erhalten wir aber bei Einführung von komplexen Größen 


Jr tineilP Pils...» : 2: 2: iR 
Ebenso wird das Impulsmoment 
L=o\rxvddr=—-olrxönds=Roe|dzdze...... 8), 


wenn R den Realteil des dahinter stehenden Ausdrucks bedeutet. Mit der Bezeichnung 
v.tir, =w, v„—ivy =m, 


U,x = i VUoy W,, Üox . iVoy Be w, 
erhalten wir ferner 


vx\ rotddr=ioml, 0X ) Bnds=Rlim,J), 


-9X\rXxrtpdr=—mz, T, —-—8u,xXUXEMNDM=-RW,2T). 
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Im ganzen ergeben sich dann für die Kraft und das Moment die Ausdrücke 


a) . u ; 
P, 7 io.J+tio wi] 0o®2,T' 
( 
(9) 
IL en : 
M, = Rem) -ioz, N] | 


in Übereinstimmung mit den von Bickley‘°) in direkter Weise abgeleiteten Formeln. 

Es läßt sich weiter zeigen, daß man die Impulsgrößen so umformen kann, daß an Stelle 
der Potentialfunktion ® das komplexe Potential 2=P+iY vorkommt. Dabei treten dann 
noch zusätzliche Größen auf, die von der Querschnittsfläche F des zylindrischen Körpers ab- 
hängen. So erhält man z. B. für den translatorischen Impuls 


J„=ioe9PQds—-owF....... ur Ba ak A 

Wenn man also ® durch die komplexe Stromfunktion (2 ersetzt, so hat man noch den Impuls 
der vom Zylinder eingeschlossen gedachten Flüssigkeit, die sich wie ein starrer Körper be- 
wegt, in Abzug zu bringen. Dabei läßt sich die Fläche F wieder durch ein Randintegral 
ermitteln, wie aus der allgemeinen Formel 

i TR os | 

\dr=— \divedr 5 Prnd 
hervorgeht, die in komplexer Schreibweise für den ebenen Fall in 


\ 
=> \(rdy yaınz)z A \ziz 


( e 
übergeht. Für den rotatorischen Impuis erhalten wir, wenn die Drehung um den Punkt 


z= 2, erfolgt 


J,=io\ ®,dz io\2Q,d2z-Zow\|z— 2,’dz 
% " C 
(ıl), 
io\ 2,dz2—iowz, F 
C 


wobei mit z, die Relativkoordinate des Schwerpunktes der Fläche F in bezug auf den Purkt z, 
als Drehpunkt bezeichnet ist. Man hat also auch in diesem Fall den Impuls der Kontur ein- 
geschlossen gedachten Flüssigkeit abzuziehen. Die entwickelte Kraftwirkung entspricht der 
zentrifugal gerichteten Trägheitswirkung dieser Flüssigkeitsmasse, die wie ein starrer Körper 
rotiert. Bei der Bildung des Momentes kann man das Potential ®, für die Drehströmung 
ohne weiteres durch die Funktion 2, ersetzen, weil das Moment der Trägheitskraft ver- 
schwindet. Man erhält dann im ganzen, wie man leicht feststellt, den Ausdruck 


L=Ro|2le)zdz+tolyvox u een ae ac 
e 


Die Anwendung der Formeln soll nun an dem Beispiel der Bewegung eines Tragflügels näher 
studiert werden. 


3. Berechnung einiger Hilfsgrößen für den Fall des Joukowskischen Tragflügels. 
Wenn wir von dem Grundkreis der {-Ebene aus- 

gehen, der durch den Punkt {= — p hindurchgeht und 7 

dessen Halbmesser a > p angenommen wird, so wird durch 

die komplexe Transformation 


ur nr DT a ne 2 


- 


‘der Kreis in ein geschlossenes Joukowskisches Profil Wach ca 
übergeführt und gleichzeitig eine konforme Abbildung der 
Außengebiete von Kreis und Profil in der z- und Z-Ebene gu ı. 
vermittelt, die für die Aufstellung der Strömung um das 
Profil verwendet werden kann (Bild 1). Als komplexe Koordinate des Kreismittelpunktes 
nehmen wir 





Abmessungen am Joukowski-Profil. 


&=beiö 


2) W.G. Bickley: Phil. Trans. A Bd. 228 (1929), S. 335; Phil. Magaz. Bd. 35 (1918), S. 500; Roy. Soc. Proc. A 
Bd. 124 (1929), S. 206; Proc. Roy. Soc. A Bd. 127 (1929), S. 186. 
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an, wobei b im allgemeinen als klein gegen a bzw. p vorausgesetzt wird, so daß gelegentlich 
die zweiten und höheren Potenzen von = bzw. > vernachlässigbar sind. Die Verbindungs- 
linie des Bildpunktes A’(£=—.p) der Hinterkante A(z2= — 2p) mit dem Mittelpunkt nennen 
wir die erste Achse, ‘die Gerade 4’0 (x-Achse) die zweite Achse des Profils. Wenn der 
Winkel zwischen beiden Achsen mit bezeichnet wird, so haben wir folgende unmittelbar 
verständlichen Beziehungen 


p+tbeidö—aei?, 
acosß=p-+tbcosö; usina=bsind,! :. * .: 2.2.2... (14). 
a@—=p’+b’+2pbcosö 


Wie in den Arbeiten von Bickley gezeigt wird, läßt sich durch eine weitere Transformation 
das Äußere des Kreises, bzw. Profils auf die obere Halbebene £’(n’ >0) abbilden. Da aber 
die Rechnungen mit dieser doppelten Abbildung einen ziemlich umständlichen und unüber- 
sichtlichen Charakter annehmen, werden wir für die folgenden Entwicklungen unmittelbar an 
die Verhältnisse in der {-Ebene anknüpfen. 

Wie bereits aus den allgemeinen Darlegungen hervorgeht, kommen für die Aufstellung 
der dynamischen Verhältnisse einige rein geometrische Größen, z. B. Flächeninhalt und 


Schwerpunktlage der Profilfläche vor, die wir zunächst ermitteln wollen. 


a) Um den Flächeninhalt der Profilfläche zu ermitteln, gehen wir aus von dem früher 
eingeführten Randintegral 
1z 
F=-—R5 Pzj2at EC BEE ee Su N 


das im vorliegenden Falle die Form einnimmt 
1 p >. 
= 2, \(e+ Fl ı-#)at. en a A A 


C=beiö+-aeir—bt,+tat, 


Setzen wir 


so wird 
‚_b, a bi+tat, ER 
PRETER it, . dl=wdl. 
Mit Benutzung der Reihenentwicklungen 
1 1 ii. PR 
ta rt ) 
1 1 bt BP 
rl Srten) 0 (16a) 
1 1 bt, Rn 
(at + bti,)? = am! - 2 re A ) 





ergibt sich 


7 


PRTETT 2. )\() 
rer) 


er 


Bi) 


m ’ i ur 2 
Wenn man den Koeffizienten der Glieder 7 mit c_, bezeichnet, so hat man nach dem 
Residuumsatz Sonn 
= =3 


Das gibt in unserem Falle für die Fläche F den Wert 
Eleelsfealt)+-] 
5 1+2(, a 17 Bu 


2. 


F=nu? 





=n.a? 


Er 
(a? — b?)? A Re a ee 
— 008 Ö 








= np? l1+(2 Nr]. — 4 f} - 
1424 sö 
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ö=> gibt den Fall des Kreisbogenprofils mit F,=0, wie zu erwarten ist. Für das sym- 


metrische Profil mit ö=0 haben wir dagegen das Maximum der Fläche bei gegebenen Werten 
von p und b 


ab 
F=4na ——. 
. (a+ b)° 
Da (a+b)" >4tab ist, so folgt allgemein 
F<nu®. 


Die Fläche F ist also jedenfalls kleiner als die Fläche des Grundkreises. 


b) Bestimmung des statischen Momentes bzw. des Schwerpunktes der 
Profilfläche. Nehmen wir als Bezugspunkt den Mittelpunkt des Grundkreises (auch Mittel- 
punkt des Profils genannt), so haben wir nach früherem 


1 f 1 { u; 
zZ 2 — 2,’ dr= > \le z.’ndf=-Zi\|z -2,%d2=(2,,—2)F. 
Nun ist mit derselben Bezeichnung wie oben 


22 =@-2)@-2)=[ pP? _ \fa, _pitt, 
z “0 =I(2 2.) (2 z)=(at+,7, b (5 ’ a E + 1). 


Benutzt man die Reihen (a) und die weitere Reihe 


\2 


1 1 bit BbY/4.\ 
EB: > 0 Re _o RR 7 
(at+bt)lat,+tbt) Kerr Bi. (>)(7 0. . (17b) 


: b ERRLE i 
und entwickelt nach Potenzen von t und — multipliziert ferner mit 


d 2__ p? u. b {A „u b | 
dal rer ) 


und integriert, so wird mit d{=adt nach dem Cauchyschen Satz 


2 e-iö 
eiö+P "5 " 


„p*b 
a? — a? — b?| 


(2,—2)f=na r 





Mit Einsetzung von F ergibt sich also 
„ »"b cos ö 
in a 
\ a@—b 
i (18). 
p® \_p*bsind 
(Ys — Yo) (1 +7 Be ”) — (a? b?)? 
Bezeichnen wir die Koordinaten in bezug auf ein in den Mittelpunkt M, verschobenes System 
mit x' und y/, so haben wir 


Be ER EN 
u hu Den an; 7 Fa zii 2,0088), 
(19). 
Be p*becosd b cos 2 3 b . ) | 
A, E ı Ba 2 (a? — 2?) (p + b cos ö) > > € » co! " 


Für ö=0 wird 
„* 
Le — RAR y,=V®. 


*72(p+2b)’ 


her T . . * ’ ® r 
Für d9=— dagegen haben wir mit «? — b’= p? die Werte 


r 


Es ist bemerkenswert, daß sich der Wert für «/, für b=0 in beiden Fällen einem von Null 
verschiedenen Grenzwert nähert, der natürlich hier nicht in Betracht kommt, da die Fläche 
des Profils verschwindet. 
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4. Aufbau der Strömung. 


Die komplexe Stromfunktion für die absolute Bewegung der Flüssigkeit infolge Trans- 
lation des Tragflügels leitet sich aus der entsprechenden Funktion für die Relativströmung 
0) 


(r)rei : ) I(C ia. ; ö ee in The use 46 


durch Hinzufügung der Parallelströmung 


ab. Man erhält dann nach Einsetzen des Wertes für z 


(21) 
So 
Einige rechnerische Schwierigkeiten macht die Herstellung der Stromfunktion, die durch die 
Rotation des Tragflügels etwa um den Mittelpunkt M, bestimmt wird. Wir haben zu dem 
Zweck zunächst das Quadrat des Abstandes eines Randpunktes des Profils vom Drehpunkt /, 
zu ermitteln. Man erhält dafür 


2 — 2,” =(2 — beid)(2 —be-iö) | 


= (a | r ‚(ae u“ iv) | 


eiy 
Tbeiö+taeit, er 


a? ,2p abcos(p +8)+2p?’a?cos2p+p* 
a?’ + b? 3a bcos(g —0Ö) 


Wenn wir unter Beachtung des Verhaltens im Unendlichen die absolute Stromfunktion (2, 
wieder in der Form 


PR © B°ePR ME UM Er 


” m 
L > >o 


ansetzen, so ist zu verlangen, daß der - der Funktion 


TR P2 io |2 20, 


auf dem Profilrand konstant wird. Wenn man (2, auf den Nenner N=a?+b?+2abcos(g 


bringt, so wird 
io [ ] 


2, =— v!4be iö+ ae in) + Be ip.N 


Bringt man nun zum Ausdruck, daß sich im Zähler des Gesamtausdrucks F (2,).., der Nenner N 


heraushebt, während der übrigbleibende Faktor gleich einer reellen Konstanten A wird, so 
erhält man drei Gleichungen für A und für die komplexen Größen A und B von folgender Forn 


bB— p’a?e-iö=V, 
Br nt 
Aa+F, e-id+Babieiö=), 
1 j R Ps. ER ; 
5 p*—be-id.A— p’a?e-?id— i(a +b)=0 


Daraus leiten sich dann folgende Werte ab für A, B und } 


ia= — / F 
2 (a? — b?)’ a’ — 


4 2? 4 
p ran "neid, | 


Man bestätigt dann leicht, daß damit auch die Geschwindigkeitsbedingung erfüllt ist. Als 
gesuchte Stromfunktion für die Absolutbewegung der Flüssigkeit infolge Drehung des Trag- 
flügels erhält man daher 

p* beid A 
(),= = er 
U EZ RSE A) | 


27). 
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Zu den beiden durch Trans- 
lJation und Rotation entstehen- 
den Strömungen tritt nun 
weiter die zirkulatorische 
Strömung um den Tragflügel 
hinzu, die der Stromfunktion 


i Be en 
pr In(£—L,) . (28) 
entspricht. Die Größe der 
Zirkulation /' bestimmen wir 
aus der Kuttaschen Be- 
dingunz des glatten Abiflie- 
ßens an der Hinterkante des 
Flügels oder des Zusammen- 
fallens des Bildpunktes der 
Profilspitze mit einem Stau- 
punkt der. Kreisströmung, 
d. h. 


>. 


Ohne Rotation erhalten wir, 
wie bekannt, den Teilwert 
I;,=4nrv,asin(a+ ß) (29). 

Bild ?2. Absolute Stromlinien (Geschwindigkeitslinien) in der Umgebung eines 
Bei Berücksichtieunge der translatorisch bewegten Tragflügels. 
Drehung wird 
(a? b’) » b(peiö-b e?id) 
| BET EN Pe PIE PER ‚hs. il ih 
h a” b? | 


l,=2ro 


oder bei Benutzung der Formeln (2) 


„, pb’cosö E b’ cos ö 
bh? = Rn h: —=4ıw IP’+ 2p b cos d u 
da“ ä | 


E- 4 20 Ja? (30). 


p+2bcos ö) 
Man sieht dann, daß damit auch die imaginären Ausdrücke in der Bedingungsgleichung ver- 
schwinden. Damit ergibt sich als Gesamtzirkulation 

‚„..pb” cos Ö\] 


"=4a\r,asin(a+P)+wla’ — bi} a (31). 


I ' 


Für die vollständige Absolutströmung in der Umgebung eines Traeflügels erhalten wir also 
> . 
die komplexe Stromfunktion 
Ba 4 p’ b eid 
io yp°| . Ari 
I ta” b’)£ 
(32). 
, j A ‚b° cos Ö\| 
I, asınla+ P)+ la? b° +1 - 2 | In (£ 
| a“ b? J 


Um eine Anschauung des Strombildes zu vermitteln, haben wir wie früher im Fall der Platte 

die absoluten Feldlinien für die Translation und Rotation getrennt graphisch dargestellt. Bild 2 

bezieht sich auf die absolute Strömung um einen Tragflügel bei Translation mit Kutta- 

Zirkulation, die mehrfach früher von mir skizziert, aber bisher noch nicht punktweise genau 

konstruiert wurde. Die Stromfunktion läßt sich nach Weglassung der zweiten und höheren 
R 


h 
Potenzen von und mit Einführung von r = 2. also 
( » ( 
Reiv=zarei" 
in die Form 
1 5; a "bi ö 
pP, -sin(y + da) 1 ; — sinig a)+ ! ip sin (29 a)—2sin(a+p)lgr (33) 
Fe - 


"„a r a’ r A 


bringen. Die Stromfunktion Y, der durch Drehung verursachten Bewegung lautet dagegen 
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on 1 
7 u la(a? — b?) 


> 


[ di: : a | 
E_- cos (p — Ö) — = ; 3 6082 (g ö) | 
a,b’a’r 


(34). 
RE 1] h 
. cos (3 q od) . r] 2 [ 1 —_ p cos Ö lg 4 | 


Die Übertragung der Punkte in den r@-Koordinaten der £-Ebene auf die z-Ebene geschieht 


mit Hilfe der Gleichungen 
x pP? \ b| p’ 1 , 
ni 3 P — [cos 00812 ))|, 
a \ ra a BEREIT | 


u . p “ b ind pP’ 1 in®o— 6 
a il pe pe sınp + „jan + er sın (29 en 
Die entsprechenden Stromlinien sind in Bild 3 dargestellt’). 
5. Berechnung der Impuls- 
größen. 


Die gegebenen Grundfor- 
meln (10) lassen sich ohne wei- 
teres auf den vorliegenden Fall 
anwenden. So hat der von der 
Translation herrührende Impuls- 
vektor die Form 


? dz 
n=i09 2,0 ar 
üs (36), 


- 0% ia. F | 


o 
wo für F der bereits ermittelte 
Wert für die Profilfläche ein- 
zusetzen ist. Man erhält nach 
dem Residuumsatz 


0P 2,0 FR 


Bild 3. Absolute Stromlinien in der Umgebung eines rotierenden Trag “ 
flügels mit Zirkulation. = P 








zov,[peit —a?e-i 


f 


Daher wird im ganzen 
J=—2nov,[p? ei — u e-ie]— ov,e-it.na’|1— — 5 
t e vlt ] A (@ + b)| 
(37) 
RER 8 
> . 2 pia , 2 ia 
—2nov,p?ei nov,ure-ie il z | 
ovV1 + A. | 77, b?)| 
in Übereinstimmung mit dem von Bickley gegebenen Ausdruck. Für die Komponenten 


erhält man 
A 


ö u a Bi BR 4 2 
J.= n0v,608a \@ + pp —2p h' 


. o [ A 0) 
J,=—-nov,Sina je l +5: TEIL Bag 2 pr | 


wo a? durch »” +b5b?+2pbcosö und a? —b? durch p?+2pbecosö ersetzt werden kann. 
Ebenso läßt sich der Beitrag der Rotation zum Impuls berechnen. Man hat naclı früherem 


dz 


u 2, „ ‚ " 
nit Fieliwile ae, na, A rn 


Für das erste Integral erhält man bei Benutzung des Ausdrucks 


} 5 beid 
io \ 3 - 2irowp?* _— (40). 
Z er FR 


3) Die Ausführung der Zeichnungen nach meinen genau berechneten Vorlagen verdanke ich Herın Kara- 


Eneff- München. 
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Das zweite Integral ist bis auf den Faktor io bereits berechnet. Man hat 
r Y a?b i p? e-id] 
J, Be ino ap a an IE + u Ge. Ran 
Im ganzen wird daher 
2bheid a? b p’e-id 
J,=-i: ’ tt 55 [ei ö 42). 
: u size ra 5) un 
Schließlich sind noch die Impulsmomente L; und L, für die Translation und Rotation 
nach derselben Integrationsmethode zu ermitteln. Wir erhalten z. B. 


i Be 
M=Re\ 2OE-2)4-dL+0uX\rdr. Er a 00 0 8 A 


Das erste Integral kann nach dem früher angegebenen Verfahren ausgewertet werden, wenn 
man etwa die Variable e#=t einführi und eid®=1t, setzt und die Partialbruchzerlegung aus- 
führt. Man erhält auf diese Weise 


p* sin (ö+a) p sin(ö+.a) , p'sin(a — Ö) 


(a? — b? a? (a? — B*) au (a? — d?)® (44). 


(L),=—-2nrov,a?b 


Der zweite Ausdruck wird bei Verwendung der oben berechneten Werte für tr, und F 


. 2 nm 
zov,a®p 
(a? — b?)’ 


(L.),.= dsin(ö+a+psinalle —b— pP?) . ” .. . (8), 
wo a?— b?—p? auch durch 2pbcosö ersetzt werden kann. Damit wird insgesamt nach 
Vereinfachung 

sin(ö+a) 2sin(ö+a) p’sin(ö—a) 


(a? — b?)? a? (a? — b?) (a? — b?)? (46) 


— . 2 4 
L:=—nov,a’bp 


Man sieht, daß L; verschwindet für b=0, d. h. für den Fall der ebenen Platte. Für ö = 


d. h. für die kreisförmig gebogene Platte ergibt sich wegen p’= a? — b’ 
Beine: » 2» 2 ra 
Der Bezugspunkt liegt in der Mitte des Mittellotes der Sehne bis zum Scheitelpunkt des 


Kreisbogens (Bild 6). Wenn man den Halbmesser des Kreisbogens mit r und seine Winkel- 
öffnung mit 2, bezeichnet, so hat man die geometrisch unmittelbar verständlichen Beziehungen 
WE NE EM u u 
a % 7 , a=rsinz, b=rsin, p=rsin; 


Daher kann man für das Impulsmoment auch schreiben 
Lean, rat zo Zoe. -. -. .» .». .» 22... MB. 
Für das durch Drehung verursachte Impulsmoment erhält man, weil das Moment der inneren 
Flüssigkeitsmenge verschwindet, 
Pf „8 i 2 02 
\ » beid a?» BE U 
L=+Reio| ar > En @—-2)33:4L...... 
. (a — b?)L % -ce) di? 


Nach einer einfachen Rechnung erhalten wir das Ergebnis 


L,= 


a? — b? 


| „ Ppb'cos2d pta?b?] 
un Pr pn 2 


u (a? —°}| 
Für b=0, ö=0 geht der Ausdruck über in den für die ebene Platte gültigen Wert 
L,=2nowa* 
der sich in Übereinstimmung befindet mit der von Omara durchgeführten Rechnung ‘). 
Für die unendlich dünne kreisförmig gebogene Platte (ö- 5) wird 
I,=2nowp* 
oder mit Einführung des Kreishalbmessers r und des Öffnungswinkels , 
.) a 


l,=2noar' sin 5 c08 


=—noomr*sin?’y 
8 . f) 


I 
> 
pn ( 


4, Vgl. M. A.Omara: Philos. Magaz. Bd. 27 (1939), S. 207. 
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6.. Bestimmung der Kräfte bei konstanter Geschwindigkeit. 


Wenn wir zunächst von der Beschleunigung ganz absehen, so ergeben sich aus den 
allgemeinen Formeln für die Kraft und das Moment in bezug auf den Profilmittelpunkt 
folgende Ausdrücke 


Pc=tovu,e ir !=4inoo)e-ir asin (a +, 

M=2aorp* sin2ß 
Mit Hilfe der Gleichungen 

A=P,mwa-+P,cse; W=P,sma—P,cosa .. .  ' 
erhalten wir daraus für den Auftrieb senkrecht zur Bewegungsrichtung und für den Wider- 
stand entgegen dieser Richtung die Formeln 

A troav;sin(a-+ P); F a 2 i ..; Ab 

mit den Beiwerten 


( r u 
2n — sın(a+P); u EZ 2 A 
pP { 


die mit den bekannten Ergebnissen der Kuttaschen Theorie übereinstimmen. Wenn man 
das Moment auf das Vorderende bezieht, das durch den Punkt mit den Koordinaten 2p,0 
dargestellt sein möge, der bei dünnen Profilen dem Kopfende sehr nahe liegt, so hat man 
für das Moment 


M'=M-2p(Acosa+ Wsin+bAcos(ö+a)+b W'sin (d +4 Er 
oder für den Absolutbetrag des entsprechenden Beiwertes 


1 lb 


Cut * Cm u 2 | 9 cos a n p cos (ö n 2 a)| Ca 


Mit Hilfe der ohne weiteres verständlichen Beziehung 
asin(a+P)=psina+bsin(d-+a) 
wird daraus mit 


zsin®Da-+tne | 


ce —- 
(Ct | 


1 


; Ei, o | ‚ 
sindö(l-+cos’a)+ { sind sin>a| { nesun2lö+ta). 
l 1 | 


4 


Wenn man bei kleinem a, cos’a »1 setzt und esin ösina, ebenso das Glied mit &? vernach- 
lässigt, so erhält man angenähert 


(ec) > 1? sin?2a+tzesind 


| 
j? (a+ pP) > nß, 


RER 
n2a+7—sınp >» 
pP 


; | 
(dt > Cmot Ca 


I 





Es ergibt sich also im Einklang mit den Ergebnissen der Messung ein angenähert gerad- 
liniger Verlauf der (e/,c.)-Kurve, deren Steigung von den Abmessungen des Profils nur un- 
wesentlich abhängt und ungefähr den Wert e„:c/,=4:1 besitzt. 


7. Einfluß der Bahnbeschleunigung auf die Kräfte am Tragflügel. 

Wenn wir zunächst von der Drehbewegung absehen und nur die tangentiale oder Bahn- 
beschleunigung berücksicktigen, so ist dem Ausdruck für die Druckkraft noch die negative 
Ableitung des Impulses hinzuzufügen. Wir haben also nach (37) 

0.J 


} > . —i0 7! 
Pc=iov,e I dt 


ee -% 


(a ] 


tiroav),e”'"sin(ia+P) 


woraus sich die Komponentenwerte 
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P,=4nrovjasin(a+Pfp)sina+rov,cosa|2p? 


P,=4nrovjasin(a+P)cosa+tror,sina|2p’—+ a’il+ 
$ _ en ta“ 


ableiten. Für den Auftrieb und den Widerstand erhalten wir also aus (54) die Formeln 


A=4nroavisinla+ß)+2rop’v,sin2a, | 


1 2 * \ ] 
W=2nop’i, | 3 il +8 £ v3) cos2a) 


Das Moment in bezug auf den Mittelpunkt des Profiles wird 


VL: 


h) ER 
M=2nrov; p?sin2a Nr 


(62). 


, . AFr . ” nd r . rY .. . ) [2 
In den weitaus meisten Fällen wird das für die Kraftwerte maßgebende Verhältnis En -g0 


klein ausfallen, daß sich die von der Bahnbeschleunigung herrührenden Bestandteile der Kräfte 
vernachlässigen lassen. Immerhin sind Bewegungsvorgänge möglich (wie etwa beim Start 
oder beim Sturzflug), wo der Einfluß der Beschleunigung ins Gewicht fallen kann, so daß eine 
grundsätzliche Beschäftigung mit den einschlägigen Fragen sich als zweckmäßig erweist. 
Wir wollen daher die durch Beschleunigung hervdrgerufenen Zusatzkräfte allein berück- 
sichtigen und entsprechende dimensionslose Beiwerte einführen. gemäß den Gleichungen 

1 


2op’v, 


W [im ’ ] 
>. ud 1+ p j* cos 2a) | 


2op’v, [2 p? (a? — b?)? 


—=nsin2a, 


R,= 


kw an; 


Die Abhängigkeit dieser Werte vom Anstellwinkel hat große Ähnlichkeit mit dem Verlauf 
der Meßwerte für die auf die konstante Bahngeschwindigkeit bezogenen Zahlen c„ und c,„. 
Während k, für kleine a etwa geradlinig ansteigt, wächst k„ vom kleinsten Wert bei a=0 
an, nach beiden Seiten, d. h. für positive und negative Werte von a in gleichem Maße. 
Schreibt man k„ in der Form 


Ku =n(y— cos2 a), 


so sieht man, daß die „Polarkurve* ein Kreis wird, dessen Mittelpunkt auf der k„-Achse in 
dem Punkt qx(>n) liegt und dessen Halbmesser die Größe z hat. Wir haben in den 
Bildern 4 und 5 die Abhängigkeit der Beiwerte vom Anstellwinkel und den praktisch brauch- 








0 VEEEEEn \ 


un 


I VEREEERE VEREIN VE 


I) 
Sı- BED 


BEN, VERREREN: 





A Due Be 





Bild 4 (oben). Beiwerte der Zusatzkräfte infolge 
Bahnbeschleunigung. 





Bild 5 (rechts). Polardiagramm für die Kräfte 
am Tragflügel infolge Bahnbeschleunigung. 
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baren Teil der Polarkurve für die Profilkonstanten a=5, b=0,707, ö =7 (p = 4,475) dar- 
gestellt. Man sieht jedenfalls, daß eine positive Bahnbeschleunigung bei normalen Anstell- 
winkeln eine Vergrößerung sowohl des Auftriebs wie des Widerstandes zur Folge hat, während 
eine Verzögerung im allgemeinen zur Verkleinerung beider Kräfte beiträgt. Bei negativem «a 
wird k„ negativ. So entsteht z. B. auch beim beschleunigten Sturzflug (c„=0, a=— ß) eine 
wenn auch kleine Seitenkraft mit 

k, = — sin 2B, 


die unter Umständen zur Änderung der Bahnrichtung beitragen kann. Interessant ist vor 
allem die Frage nach der Lage der Kraft, bzw. der Verschiebung der Kraftlage infolge der 
Beschleunigung. Wir benutzen den Ausdruck für das Impulsmoment, aus dem sich das Moment 


W- odL, _rot,a’bp*/sin(ö+a) 2sin(ö+a) p*sin(ö— a) 


ot a? — b? a?’ — b’ a? (a? — 5’ we 


| 


ergibt. Daraus leiten wir die Gleichung der Wirkungslinie der nur durch die Beschlennigung ;, 
verursachten Kraft ab. Mit den Werten (60) erhalten wir mit v,—=0 


9m? __ Mn | | p' 
2p a (Ita PIE 


a?bp* sin(ö+a) 2sin(ö+a) p’sin(ö—.a) | 


f 


u. 
zsina|2p’+ ax (1 ta E 5) ycosa 


a—b?| a —#? a? (a? — b?)? 


Differenzieren wir nach a, so haben wir zur Bestimmung des Schnittpunktes zweier benach- 
barter Angriffslinien die Gleichungen 


en Vai a’byp* ee (ö+a) p sin(ö-+ a) | 


a—b| @® — a? (a? — 5? 


a? b p' E (ö+a) 2ecos(ö+a) p? cos (d — a) | 


xcosaA+ysinaB=—— —— — - = Ar 
rYy a —bi a —W a? (a? — b? 


wo zur Abkürzung 


gesetzt wurde. Eliminiert man a, so kommt 


a?’byp* 1 { p* 
(= 5 hucoso| = _ 12 ; + 2  b2j? 


7} = l, 
a —b a? j? (a 


2} ’’ 1 ‘ 2 ’ 
yB= —E sin ö PPC = + ? | —=sindöD | 


a_b Ar (a? — b?)? 


Es ergibt sich also das bemerkenswerte Ergebnis, daß die Wirkungslinien der Kräfte für alle 
Anstellwinkel durch denselben invariablen Punkt hindurchgehen, der in bezug auf den Mittel- 
punkt des Profils die Koordinaten 
D sin ö 
Y ar B 
hat. Wir wollen diesen Punkt kurz als Beschleunigungspunkt bezeichnen. 


Betrachten wir zur Prüfung dieses Ergebnisses etwa den ausgearteten Fall des unend- 
lich dünnen kreisbogenförmigen Profils, so ergibt die Rechnung bei Benutzung der Beziehungen 


(65) 


a |—-—--cos? a) = #1 —- 02a 
>” / er } 
1 cos’ 


ame Be 


Die Koordinaten des Beschleunigungspunktes werden 


7” 
5 . 


Der Abstand des Punktes vom Kreisbogen wird daher 
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) dar- 


ıstell- 
hrend 
em a 
) eine 


t vor 
e der 
ment 


(64) 


Ing ı, 


nach- 


alle 
ittel- 


(65) 


end- 
gen 
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Der Beschleunigungspunkt fällt also im Falle des 
Kreisbogens mit dem Kreismittelpunkt zusammen?). 
Die durch diesen Punkt hindurchgehenden Kräfte stehen also 
senkrecht auf dem Kreis (Bild 6). Die Gleichung der Angriffs- 
linie erhält die Form 
xsinß(p+a’)+ycosab’+p’beosa=0 . . . (69). 
Transformiert man auf den Kreismittelpunkt, so erhält man 
die reduzierte Gleichung 
zsina(p+a)+ycosab®’=0. . . . . . (69a). 


u % 
Bild 6. Zur beschleunigten Translation einer kreisförmig gekrümmten Platte. un 4 


8. Einfluß der Drehung um die Querachse auf die Kräfte am Tragflügel. 

.Wir wollen jetzt von der Bahnbeschleunigung absehen und eine gleichmäßige Dreh- 
bewegung um die Querachse voraussetzen, die der konstanten Winkelgeschwindigkeit ® ent- 
spricht. Unsere allgemeinen Formeln geben dann ohne weiteres als Ausdruck für die auf 
den Tragflügel wirkende Druckkraft 

P=—-ioJ+tiov,e-i«l 
re. Ri LT N a? / :d p? e-iö\ 
Daraus leiten sich folgende EEE ab 
’ : b’ cos ö 

P,=4nov,sina |v,asin(a+P)+ w (a*- -b? +: . a ) 


| p* } f b ) 
-zov,wsina 29% (1 +a@- Er) - id 0 „? a le 


b° 
P,=4r0ov,cosa I" asin(a+P)+ o(@ +; Ä ee) 


| ) „nee ptin 24 “ ( x na 
— b?}? —b a? — b? a? — b? 


a? 


7„* 





+rzov,wcosa|2p? -e(14 


die dann weiter Auftrieb und Widerstand ergeben in der Form 


j b’ cosö a 
A=4novsasin(a +) +rov,w 1{@- nr ) 2 p’cos 2a— 014, 7) 


2 24 ’ 
> r (3a?b — 2b’)sin(ö+ a) — Be sin(ö—a)|, 


zoa’p' 


W=2rov,wp’sin2a+ — de (3a?b — 2b?) cos(Öö+«a)+ p 





a? = E cos (d — a) 
Das Moment der Kräfte in bezug auf den Profilmittelpunkt wird ferner gemäß den allgemeinen 
Grundgesetzen 

M=R(iv,ei“.J) 


2 ms 
—2novp’sin2a+rov, ut — b?? 30 b—23°)cos(dö+4)+-; cos (0 a| ö 


Setzen wir wieder 


1 1 1 
A=ta zom ip; Wen goun dp; M=emzeonltp), 


so entstehen schließlich * die dimensionslosen Beiwerte die Formeln 
a opf,[®—b , beosö] „ = p 
.=2n „ sin(a+P)+ z* r 4° p’ ur, (a — —_ +2cos2a — 1+- (a’ en) 
1 »p 2 12 
fr 2, ) (a? - Bar 


(3a? b — 2°) sin (6 +0) sin (d-a)l, 


= Psinta+ga(®?) zur (3a?b — 2b?) cos(ö + a)+, 5 cos (d—a) ; 


2 — b?)? 


®, 


2 n? b 
u cos (ö — a) 





8 


5), Vgl. W. G. Bicekley: Philos. Magaz. a. Journ. of Se. Bd. 55, 6. Ser. (1918), S. 396; Philos. Trans., Ser. A, 
Bd. 228 (1929), S. 2335. 


m=—asin2a+4n GE p ji Bad — 28) cos(d+a) + 
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b 
Wenn man : e einführt, so wird z. B. angenähert bei dünnen Profilen 
) 


p' p? i i 
2 u VECOSÖ), 
a? b? ° 1+2ecosö pl E COS ( 
3b — 21° 3(»2?-+2pbeos od) + b? 3(1+28cosö)-+ e? 
p' - 3.» [2 p' 1 Pr f > u € p’ ° — ER 3: pi Dr cos d). 
Br (p?+2pb.cos ö)’ (1 +2* cos Ö)? 
a? bh ep’ (1l+28cos Ö) 
0 0) + rs 3 - S ) 
N (a* b?)? (1 -+2ecos Ö) ep (1 te cos 0) 


Bei Beschränkung auf die erste Potenz von # kommt daher 


en Ad np z op \ { 
a =2n r sin (a 4) +2a| | (eos? a+2#c0s od) a| E [2cosdsina + sindcosa], 
) ’ v® 
0 0 
i 2 wp u in) p “ RR N er r s r ER 
„=a|- |sin2a+n s | e[2 cos öcos a — sin Ösin a], (75). 
0 \ 0 
TEE DER, 9 
mM=y7 sin2za ga e [3 cos (dö+a)+cos(ö - a)] | 


Wenn man das Moment auf den Punkt 2p,0 bezieht, so hat man nach der unter (56) auf- 
gestellten Reduktionsformel nach gleichzeitiger Änderung des Vorzeichens für den Absolut- 
betrag des Beiwertes 





= (le } 
E* : m en 5 n a 8,5668 er 4 | eis 4% en +2 - pP a)\ 
= n ve cosa— 2r u > cosd —ı ir a ja) COS Öcosu, (76). 
+ n ve . coS a | 1 + a: ea) 
„En : ; (”P\sin Ö - ID b En 1 un - cos o) 





wo der nur von der Translation abhängige Teil (c/,); bereits früher näher bestimmt wurde. 
ER j B b 
Bei Vernachlässigung der zweiten und höheren Potenzen von r] —e erhält man daraus 
. E non a 1 wp * . R 1 op\” } 
C„=ld, ht — 7 P cosa - y?' P; [4 cosöcosa + sinösina] +5 a| ‚) esind. 
’o ®, ”, 
Zur Veranschaulichung dieser Formeln haben wir die Beiwerte für einige Werte des 
.- u ‚’ pP. u . . . > R e £ 
Verhältnisses a Abhängigkeit von dem Anstellwinkel aufgezeichnet. Im gewöhnlichen 


0 
Bereich der Anstellwinkel verlaufen die c,- und c„-Kurven ungefähr geradlinig. Während 














N 2 P r . . . “ ap 
aber die Auftriebskurven angenähert parallel liegen bei den verschiedenen Werten von , 
0 
12-—— a0 — ———— 11 
id— 016} N N N 1 
1 | 
’ 
2 —+ ne GN BB rn [up re 
| Bild 7 (links). | 
27 Einfinß der m — rt 
| Drehung eines 
Tragflügels 004 + mu + + 


| auf die Auf- 
TT71 triebswerte. ” 05 47 
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24 
A ae a Er DR SER EEE 
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— Bild 8. Einfluß der Drehung eines Tragflügels auf die Wider 


standswerte. 
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derart, daß c„ wächst bei positiver Drehung, so erfahren die Widerstandskurven bei An- 
wachsen des Verhältnisses eine Schwenkung im positiven Sinne oder eine Vergrößerung der 
Steigung 


7. Einfluß der Drehbeschleunigung auf die Kräfte am Tragflügel. 


Wenn wir bei der Drehung eine Tangentialbeschleunigung annehmen, so sind die da- 
durch hervorgerufenen Kräfte durch die Beziehung 


3J, [2 beid 25  e-iö 
OJr . , IS0e da , pe 
P ae =inowap'|— ‚tr » ;[eid 7 =" 77 
of ee. la?—b? (a? —b?) a? )\ A 
bestimmt, aus der sich dann wie früher folgende Auftriebs- und Widerstandswerte 
' aur 13 a?b 2b° s a?b p” 5 | 
A=no0owp| 5 3,5 cos (c at 55 COS (( a)|, 
| (a? — b?) \ (a? — b?)? | 
Rn (78) 
3a? b— 2b‘ a?bp* | 
W=rzowp* | = =— sin (Ö a) = >; sin (Ö a)| 
a 1 L ka’ b?)? s + (a’ - b?)? | 
ableiten, denen entsprechende Beiwerte, etwa 
A W 
a 77 “Toop* 
zugeordnet werden können. Da das Moment durch den Ausdruck 
Y OL, 2rowp*| „ p’b*cos2ö p* a? b? ] "0 
M=- ” r a 1 — 7} 2 = Fr J) 
Of a—b: | (a? — b??) (a? — b?P) X 


bestimmt ist, so erhält man z. B. die Angriffslinie der durch Beschleunigung hervorgerufenen 
Zusatzkraft, in der Form 


RE 08: p? | An u. o* \ 
ee lt p)|+ysind at gli >) | E 


(SO). 


a u 


p’b'cos2ö6 pra?b? | 
(a? - b?)? (a? . b?°| 





Wir wollen in diesem Zusammenhange auf weitere Einzelheiten verzichten und nur noch den 


Spezialfall der kreisförmig gebogenen Platte FR. p=a’— b°) kurz erwähnen, bei dem 


die Impulsgrößen die Werte 


@% 
J,=2nowp?b; L,=2nrowp* a 3 + r 
annehmen. Da ferner die letzte Gleichung \ j 
\ 
2 a 
1. N 
liefert, d. h. die Ordinate des Mittelpunktes des Kreisbogens in * 
bezug auf das in den Profilmittelpunkt als Drehpunkt verlegte m \ly N 


Achsensystem, so ergibt sich, daß die Kraft durch das Zentrum 
r . . = r . i 9, 7, Se y ipte E 
des Kreisbogens hindurchgeht und der Sehne des Kreises parallel Pld 9. Zur beschleunigten Ro 
B 5 . tation einer kreisförmig ge- 
gerichtet 1st. krümmten Platte. 


8. Kleine Drehschwingungen des Tragflügels um die Querachse. 
Um schließlich die Wirkung kleiner Drehschwingungen des Flügels um eine mittlere 
Lage des Anstellwinkels zu untersuchen, setzen wir 


o=w,sin?n T 


und bestimmen die über die Zeitperiode T der Schwingungen genommenen Mittelwerte aller 
Ausdrücke, welche die Geschwindigkeit und Beschleunigung enthalten. Dabei geben nur die 
quadratischen Glieder von Null verschiedene Beträge. Wegen 
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ergeben sich dann folgende Beiwerte für Auftrieb und Widerstand 


3a?b— 2b’ vb . 
den r 5 — u sin(ö — a) 


> 
- (a? — b?)? 


G> sin(ö+ «) - 


e 
n 3n ein (a+ Pß) 


D.P\?[_: 2 var | 
7 P\ p nr cos (+) + pi; cos (d —a)| 


vo.) - b?)? 


/ 2 
[73 ) ö; . 
| R\ € (2 cos a cosdö —sina sind), 
o 


\2 


|-esin 6. 





@®,Pp 
x of 


, 
\ % 


Man sieht, daß kleine Schwingungen im allgemeinen eine Verkleinerung des Auftriebs und 
eine Vergrößerung des Widerstandes zur Folge haben. Wenn man bedenkt, daß Schwingungen 
dieser Art als Begleiterscheinung der periodisch sich abspielenden Wirbelablösung auf der 
Hinterseite des Flügels vorkommen, so ist damit zum mindesten eine gewisse Seite des 
Widerstandes theoretisch verständlich gemacht oder eine theoretische Bestätigung dafür ge- 
geben, daß der Formwiderstand durch die oszillatorische Beschaffenheit der Strömung nicht 
unwesentlich mitbedingt ist. Eine Unstimmigkeit gegenüber der Erfahrung besteht allerdings 
darin, daß die Widerstandszahl ec, unter den gegebenen Voraussetzungen mit dem Anstell- 
winkel nur wenig veränderlich ist und mit wachsendem «a sogar etwas abnimmt. Man muß 
aber dabei bedenken, daß die Kräfte auch noch durch die Lage des Drehpunktes bestimmt 
werden, der in unserem Falle in den Profilmittelpunkt verlegt wurde. Eine weitergehende 
Verfolgung der Theorie müßte auch diese Abhängigkeit in Rechnung stellen, um eine weitere 
Annäherung an die wirklichen Verhältnisse zu erreichen. Als physikalisch maßgebende Größe 
geht ferner das Verhältnis der Drehgeschwindigkeit zur Bewegungsgeschwindigkeit in die 
Formeln ein, das natürlich empirisch ermittelt werden müßte, um eine Prüfung der abgeleiteten 
Ergebnisse zu ermöglichen. Es scheint, daß diese Schwingungsfrequenz für die vollständige 
Erfassung der Kräfte am Tragflügel nicht entbehrt werden kann. 


| TTTTTT7 Der besseren Übersicht halber haben wir noch das 
JEREHEBERSERE Polardiagramm für einen Joukowski-schen Tragflügel 





’ r n a a 
mit den Abmessungen a=5, b=0,71, ö =T aufgezeichnet, 


Bild 10. Dabei wurde das Geschwindigkeitsverhältnis 





t 2 
(®8 P) =0,05 zugrunde gelegt und außerdem noch der 
\ 0 

Einfluß der endlichen Länge des Flügels, d. h. der induzierte 
Widerstand 

















2 
Ca 
Cwi= A 











und die dadurch hervorgerufene Verkleinerung des c„- Wertes 
berücksichtigt, der für die Flügelstreckung A bei Vernach- 
lässigung des Einflusses der Drehung den Wert 























2rsin(@a+P) 
. F} * .. Ca Fr > e 
Bild 10. Polardiagramm eines Tragflügels i+ A cos(a+ P) 


N 2 
mit 46, (>?) =0,05. 
vo 


annimmt°®). Auch die mit der Messung gut übereinstimmende Abhängigkeit der Momenten- 

zahl c„ von der Auftriebszahl c, ist dem Polardiagramm hinzugefügt worden. Wegen der 

ut ist natürlich nur ein qualitativer Vergleich 
o 
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Unbestimmtheit, namentlich des Quotienten 


mit den Daten der Messung zulässig. 


6) Vgl. Wilhelm Müller: Mathematische Strömungslehre, S. 166. Berlin 1928. 
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Die Grundgleichungen der elastischen Stabilität 
in allgemeinen Koordinaten und ihre Integration'). 


Von H. Neuber in Braunschweig. 


Ausgehend vom Kräftegleichgewicht stellt Verfasser unter Anwendung des 
absoluten Differentialkalküls die Grundgleichungen der Stabilitätstheorie 
in beliebigen krummlinigen Koordinaten auf. Der Vorteil dieser Ableitung 
ist der, daß zunächst keine Voraussetzungen über das Formänderungsgesetz 
gemacht zu werden brauchen und daf die Größen, denen auch bei end- 
lichen Verschiebungen Tensoreigenschaft zukommt, klar als Tensoren her- 
vortreten. Das Integral der Grundgleichungen des anisotropen Kontinuums 
bei stationärer Vorspannung wird auf eine vektorielle Spannungsfunktion 
zurückgeführt, die einer Differentialgleichung sechster Ordnung genügt. Im 
Falle der Isotropie gelingt die Zurückführung auf Differentialgleichungen 
zweiter Ordnung. 


1. Einführung. 

In den seitherigen Darstellungen der Stabilitätstheorie [1 bis 8]?) werden zur Aufstellung 
der Grundgleichungen kartesische Koordinaten verwandt. In vorliegender Arbeit wird. diese 
Theorie durch Anwendung des absoluten Differentialkalküls in beliebigen krummlinigen Ko- 
ordinaten hergeleitet. Die Darstellung bietet den Vorteil, daß jene Größen, denen auch bei 
endlichen Verschiebungen noch Tensoreigenschaft zukommt, besonders klar als Tensoren her- 
vortreten. Ein wesentlicher Unterschied gegenüber früheren Darstellungen besteht vor allem 
darin, daß die Gleichgewichtsbedingungen nicht aus dem Energieprinzip, sondern allein aus 
dem Kräftegleichgewicht hergeleitet werden, d. h. noch ohne eine Voraussetzung über das 
Formänderungsgesetz bzw. das elastische Potential. Nach Elimination des zusätzlichen, durch 
die elastische Deformation bedingten Spannungstensors mit Hilfe der Spannungs- Dehnungs- 
gleichungen ergibt sich für den Verschiebungsvektor der anisotropen Kontinua ein neues 
Gleichungssystem, welches einer Modifikation des anisotropen Elastizitätsgesetzes entspricht, 
derart, daß an Stelle des ursprünglichen Elastizitätstensors ein neuer Tensor tritt, der von 
den Vorspannungen abhängt. Das allgemeine Integral der Grundgleichungen bei stationären 
Vorspannungen wird auf eine vektorielle Spannungsfunktion zurückgeführt, welche einer 
Differentialgleichung sechster Ordnung genügt. Im Falle des isotropen Kontinuums läßt sich 
zeigen, wie der gleiche Ansatz, der vom Verfasser in der gewöhnlichen räumlichen Elastizitäts- 
theorie mit Erfolg angewandt wurde, zu einem einfachen Lösungsweg führt. Hierbei treten 
Differentialgleichungen von nur zweiter Ordnung auf. 


2. Bezeichnungen und geometrische Grundbeziehungen. 

Es mögen mit z=x', y=x°, 2=x°, allgemein x* orthogonale kartesische Koordinaten 
gekennzeichnet sein, und es mögen römische Indizes stets auf kartesische Koordinaten hin- 
weisen. Die Größen x* sind als differenzierbare Funktionen der neuen krummlinigen Ko- 
ordinaten x', x?, x’, allgemein x* vorausgesetzt, und es mögen griechische Indizes stets auf 
diese krummlinigen Koordinaten hinweisen. Jeder, die Komponenten eines Vektors in karte- 
sischen Koordinaten darstellenden Größe Ak= A, ist in den neuen Koordinaten die Größe 


3: 
A = > Ax dr! . . . . R “ s ö . > (1) 
k 


zugeordnet, die auch kurz als Vektor oder Tensor erster Stufe bezeichnet werden kann. Diese 
Art der Unterscheidung der kartesischen Ausgangskoordinaten von den allgemeinen Koordi- 
naten weicht zwar von der üblichen Schreibweise ab, hat aber bei der Differentiation gewisse 
Vorteile, insofern als kartesische Vektoren für die neuen Koordinaten wie Invarianten be- 
handelt werden. Im übrigen gelten die bekannten Bezeichnungen (vgl. z. B. Levi-Civita [9)). 
Wird zur Abkürzung 

our ck (2) 

ea er an aan 


geschrieben, so bilden allgemein die Größen Bi., die sich entsprechend der Beziehung 


1) Vortrag auf der Tagung der Mathematiker-Vereinigung in Würzburg 1943. 

2, Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. Während 
der Drucklegung erhielt ich Kenntnis von einer Arbeit von E. A. Deuker (Deutsche Math. Bd. 5 (1940), S. 546 bis 567), 
in welcher bereits allgemeine Koordinaten verwandt werden. Das Problem wird dort durch Einführung eines Hilfs 
parameters als Variationsaufgabe dargestellt. Die Variation ist jedoch nicht allgemein ausgeführt, sodaß die Grund- 
gleicehungen nicht explizit erscheinen. Auch ist die Integration dort noch nicht behandelt. 


+) 
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’ kl a 
Biu = #) Bxı C, Cu . . . . . . . . . . . . . (3) 
kl 
aus den Größen B;,; ergeben, einen „Tensor“ zweiter Stufe. Entsprechendes gilt für Tensoren 
höherer Stufe. Der für die Längenbestimmung maßgebliche metrische Tensor wird 


u N chel RE TE FE 
k 


; j m k : en 
mit den kartesischen Komponenten 9,,;=d,. Hierbei ist 
l, wenn k=] 
‘ 0, wenn kl 
das Kroneckersche Symbol. Bildet man ferner in der Determinante der c{ die Unter- 


determinante von ec‘ und dividiert sie durch die Determinante |c}, so ergibt sich die Größe 


BI, 1 B* 109. „m > 
= k e, EkIm® a Ba ur na wur : 
2c; 
Al l,ın,a,? 
Für die hierbei auftretenden Hilfsgrößen e;;„ gilt 

l, wenn k,l,m eine zyklische Form bilden; 

1 bei ungerader Permutation 
Hin der Indizes gegenüber oben; } 

0, wenn mindestens zwei Indizes zu- 
sammenfallen. h 


Das Entsprechende gilt für die Größen e‘“’. Für die Determinante folgt 


k N 1 kl 
: „iR ‚m ) 
C; = G kım e@ 37 cn czc, . . r r . . . ° ‘a, . (5). 


k,l.m,a, ?,% 


Nach der Determinantentlieorie gilt, wenn nunmehr stets bei Summation über einen Index, 
der einmal oben und einmal unten steht, das Summenzeichen weggelassen wird: 

»B__ 58 k„.a__ sk rn 

u #0 PU 507 SE” MB | 
Setzt man ferner 


1 
Ja 3| —— 6 e«?Y er! Je} 3“ Iyr = 4 . . . . . . . . . (10) 


und 


a Er a 


ge? - e«y! 
so folgt j 
BE 5: en Dina 


d. h. auch die g“* bilden einen Tensor. Ferner wird 
k f 1 
= Vg, c'= SEN ER HE Fra en Sr 
ER v9 
und die Größen c;, ö;, sowie 
1 
—e 
v9 
stellen sämtlich Tensoren dar (hierbei muß allerdings verabredet werden, daß das neue Ko- 
ordinatensystem in demselben Sinn dreht wie das kartesische). 
Die Größen Bi, entsprechend Gl. (3) sind kovariante Tensorkomponenten, während die 
zugeordneten Größen 


Met Bi Ygaymtlı - .» 2 2. 222. 


a WE 
die kontravarianten Komponenten angeben. Es folgt 

a u ut FRE Su 5 
Bei Differentiation einer Invariante c entsteht ein Vektor 


dc 


za aeiale : 2... vn. 
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Bei Differentiation eines Vektors A, entsteht der Tensor 


d Ar ) 


3 u. [„& 
Auszac hne 4: 
Ar = zu cr A 323% Ac) 


Die entsprechenden krummlinigen Komponenten dieses Tensors sind mithin 


ie . 
Ax,ı a dr’ I a? A; 


Hierbei ist 
k OCh „dc 


Ca ‚=z— (Cr 
0x " 


=-—T! 
0x’ P 
gesetzt (entsprechend dem Christoffel-Symbol zweiter Art). 
Allgemein folgt als Regel der kovarianten Differentiation, daß die Operationen 
Ö Ba 3 i 
Br 
a» Ö b® ‚ va ‚ 
B? = — +19, Br +09, Ber 


0x7 


Ba 3,, = r% Y Bi Au Tr; y Baı 


usw. stets neue Tensoren liefern; ihre kontravarianten Komponenten sind 
Baa|!=g'*Bas,,, BePlr= BR ,g’* usw. 


Weiter folgen hieraus die Ricei- Beziehungen 


k X aß e ö 4 
Ca 3,3 May = I... g° |? = a3]! = 0 


3. Spannungstensor und Gleichgewichtsbedingungen. 
Werden die Normalspannungen o,, o,, 0, mit #£,,, #33, #33, die Schubspannungen r,,, 
T;x; Txy mit t,,, t,,, t,, gekennzeichnet, so repräsentiert allgemein f;, die Komponenten des 
Spannungstensors in kartesischen Koordinaten. Die zugehörigen kovarianten Komponenten 
im krummlinigen System sind dann 

k 1 . 
A — txı C, Cu . . . . . . . . r . n . ’ . (25) 

bzw. die kontravarianten Komponenten 
a ee a 


Die Gleichgewichtsbedingungen gegenüber einer Verschiebung in x- Richtung 


00,  Oftxy | Org; > u 9" 
'y 7 Da T Bl ag a 7 Er ee (27) 
gehen in die Form 


ee re 


über. Hierbei ist P! der Vektor der Massenkraft. Entsprechend ergeben sich als Gleich- 
gewichtsbedingungen in den neuen Koordinaten 


er 2 MR BG EEE NE: . (9), 
wobei P'" aus 
a 1 2 


zu bestimmen ist. Die Gl. (29) lauten ausführlich 

Ri) ru 

dx 
Eine weitere Bedingung ist das Gleichgewicht der Spannungen gegenüber einer Drehung. In 
kartesischen Koordinaten folgt hieraus die bekannte Symmetrieeigenschaft k!=1t!k des 
Spannungstensors, die infolgedessen auch für #“ gilt. Zur Darstellung der physikalischen 
Spannungskomponenten 7X“, die nicht mit den Größen 1X“ identisch sind, muß der an einer 
Fläche x* = konst. wirkende Spannungsvektor ermittelt werden. Denkt man sich diese Fläche 
als vierte Seitenfläche eines Elementartetraeders, dessen übrige Seiten durch Flächen x* = konst. 
ran werden, so gilt die Gleichgewichtsbedingung für die resultierende Kraft in der &- 


ichtung 

l i 

RE - ‘ 

Wr IR. ER R. De a a TE 1) 5 
u Yg«. ; Ih 


+ FFIR PP PEN : . 2, 2: 
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Hierbei stellt 1%“ die #- Komponente des an der Fläche x* = konst. angreifenden physikalischen 
Spannungsvektors dar. c./vg« „ist Richtungskosinus der «-Richtung gegenüber der i- 
Richtung, während a/vg* Richtungskosinus der Normalen zur Fläche x* = konst. gegenüber 
der x*-Richtung darstellt. Durch Multiplikation mit c;, und Summation über 1 folgt 


a RER EEE 


e g*’ 

Die Größen ı*’ sind mithin kein Tensor. Die Normalspannung wird 
gar 

g* 2 

während die gesamte, in die «-Richtung fallende Komponente der Schubspannung die Größe 


(34), 


OL = 


a EARTH 
el VI" Geo 


Tiu 


annimmt. Sind längs der Oberfläche des Körpers die auf die Flächeneinheit bezogenen 
Kräfte f! vorgegeben, so lauten die Randbedingungen, wenn die Oberfläche durch x* = konst. 
dargestellt wird: 
ER op 
2 uymm—l. Ei ee 
4. Geometrie des deformierten Körpers. 
Nach der Deformation gehen die Koordinaten »* über in 
ee ee a 
Die Größen V* stellen die kartesischen Komponenten des elastischen Verschiebungsvektors 
dar. Werden entsprechend alle neuen Größen durch einen Strich gekennzeichnet, so folgen 


die Beziehungen 
k_kıvk_Xkı kr on 
FREE RR TE LE EEE u iR 


Yan — Giu +: V. n 2 Vu,a — Be V,,i . . . . . . . . . . (39), 
Vg=yYgil+ V’,„+ Glieder höherer Ordnung). . . . . (40). 


Ebenso lassen sich alle weiteren Größen auf die entsprechenden Größen des Ausgangssystems 
und die Ableitungen des Verschiebungsvektors zurückführen. Der Verschiebungsvektor V; 
und seine Ableitungen stellen im Rahmen der nachstehend behandelten Problemstellung für 
das elastische Kontinuum bereits Größen zweiter Ordnung dar, da die Komponenten des Vor- 
spannungstensors T*“ als Größen erster Ordnung anzusehen sind. Bei Beschränkung auf 
Glieder zweiter Ordnung ergibt sich 


Yan = (iu n 23. 7 7 Pi A LANE OR Ve 
Hierbei ist 
di u — { Vı, “ - Vi, 2) . . . . . . . . . . . (42) 


der Verzerrungstensor. Weiter erhält man 
r=4-V.,c, are... ie a 
En 2, 
Die zweiten Ableitungen des Verschiebungsvektors lassen sich entsprechend der Beziehung 
ee N Te A 


die aus Gl. (42) hervorgeht, stets auf die ersten Ableitungen des Verzerrungstensors zurückführen. 
Wird die kovariante Ableitung im deformierten System mit einem * gekennzeichnet, so folgt 


CGa=(,a, As 3 = Aa, > Vias A- N Win a (46), 


A" Bu v: 3 2 Er EN ER a (47). 
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5. Gleichgewicht des deformierten Körpers. 


Ist K*“ der Spannungstensor des deformierten Körpers, d. h. 
en , Am 
IRRE a ee a ME 


und ist P“+p“ der Vektor der Massenkraft nach der Deformation, so folgen als Gleich- 
gewichtsbedingungen aus dem Vorhergehenden außer der Symmetrieeigenschaft K*“—= K4 
die Beziehungen 

ENTER u ar a er 
bzw. 


Kr, +Vi.,Kre+V",.,Ku+Prtpe=0 ...... 0. (0). 


Die so gewonnenen Gleichungen bilden die Grundlage zur Behandlung des Stabilitätsproblems. 
Waren vor Eintreten der Verschiebungen V* bereits Vorspannungen T*“ vorhanden, so setzt 
sich der Spannungstensor K*“ aus-einem Tensor 7*“, der allein mit den Vorspannungen in 
Zusammenhang steht (Größen erster Ordnung), und einem durch die Deformation bedingten 
Tensor f“ (Größen zweiter Ordnung) zusammen: 


Kltu—_ pur Tin ee a ren rien Dr (51). 


T erbei stellt sowohl 7?“ als auch f“ einen symmetrischen Tensor dar. Zur Feststellung 
aes Zusammenhanges von T*“ mit dem Vorspannungstensor T*“ muß auf den physikalischen 
Vorgang Bezug genommen werden’). Die an einer Fläche x' = konst: tatsächlich wirkende 
Kraft in k-Richtung ist vor der Deformation 


N T!«Vgekdx’dx’ (Flächeninhalt Yyg" dx? dx’) 


[7 


und nach der Deformation 

zpr; gehkdader. 
Die projektive Operation, die durch die Multiplikation mit ch gegeben ist, bedeutet zwar die 
Bildung der exakten k-Komponente; in Wirklichkeit hat aber das Material .in der näheren 
Umgebung des betrachteten Punktes insgesamt, d. h. von der Verzerrung abgesehen, eine 
Drehung vollführt (d. h. auch als starrer Körper), die durch den schiefsymmetrischen Tensor 


1 
Mu (Vic — Vu, i) 


gegeben ist. Die Ableitungen dieses Tensors lassen sich übrigens restlos auf die Ableitungen 
des Verzerrungstensors zurückführen; es folgt nämlich mit Bezug auf Gl. (45): 


Viu,v = di Te d.. nd A us » . ern . (53). 


Werden die oben gekennzeichneten Kraftgrößen nicht unmittelbar als einander gleich an- 
gesehen, sondern erst nach Elimination der in c’* enthaltenen Dreliung, so ist an Stelle von 


m ch + ck (di, +0" „) der Tensor 
+k ; k k r 
c ee 0: A 


der nur noch eine reine Verzerrung enthält, zu setzen. Dieser Gedankengang führt, wenn zu- 


r [2 . .. ” ” ” 
nächst ein Tensor 7?“ an Stelle von 7** eingeführt wird, zu einer Beziehung von der Form 


T}u V q ck — Tr u V g' (ck En ch d,,) ‚ 
bzw. a a VE 
Tr — Tar(1+di)+ Trrdi 


3) Die exakte Berechnung des Tensors 7*'‘- würde eine Bezugnahme auf die mit dem Vorspaunungstensor durch 
das Forimäuderungsgesetz verknüpfte Vordeformation erforderlich machen. Eine solche Rechnung würde eine Erwei- 
terung der Ausgangsgleichungen auf Größen dritter Ordnung erforderlich machen, da der metrische Tensor des Aus- 
gangszustandes Größen erster Ordnung darstellt und jener der Vordeformation bereits Größen zweiter Ordnung 
enthalten würde. Zugleich wäre eine Erweiterung des Elastizitätsgesetzes auf Größen zweiter Ordnung erforderlich, 
wobei entsprechende Aussagen kanm durch experimentelle Befunde sichergestellt werden können. Bei Umgehung dieser 
Schwierigkeiten muß über das Verhalten des Vorspannungstensors bei der Deformation eine besondere Beziehung 
eingeschaltet werden, für die ich nachstehend in Anlehnung an das physikalisch zu erwartende Verhalten zwei 
Formulierungen angebe. 
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Die Auflösung nach Ti“ liefert im Rahmen der Theorie zweiter Ordnung 
EB FE . .. 2... ec 


m. . . . . . m . .n ° 
Man erkennt, daß der Tensor 7T*“ unsymmetrisch ist, also nicht direkt mit 7*“ identifiziert 
werden kann. Setzt man seinen rein symmetrischen Anteil gleich 7*”, so ergibt sich 


T+ = Ti«(1 —d)) (Tr ai + Ton 0 SR ARE 
Gl. (51) liefert nunmehr 
KH — te + Tel —d,)— 5 (DT d,+T'“d}) 
und die Gleichgewichtsbedingungen (50) gehen über in 
+ TN— n (T"d+Ted)ı+V\.T®—q4T+P"r+p"=0.. . (9). 


Hierbei ist davon Gebrauch gemacht, daß t*“ von derselben Ordnung klein ist wie V?|«. 
Wird ferner beachtet, daß sich der Vorspannungstensor 7*“ mit der ursprünglichen Massen- 
kraft P“ im Gleichgewicht befindet, d. h. 


Ber ae, (60) 
erfüllt ist, so läßt sich Gl. (60) auf die Form 
iu m; rt 1 mi 7 mn 4 D 
rn rT7V.,+-z(7t& +7 a, trP=0. . .... 


bringen. Werden alle Formänderungsgrößen durch die Ableitungen des Verschiebungsvektors 
ausgedrückt, so ergibt sich folgende Gleichgewichtsbedingung: 


Au 2 rt 1 ru r 2 rk rı3\ ’ 
tr PH V te — TI, tu) THU, +9) Tr) + Pra=0 a 3 


Werden andererseits die Ableitungen des Verschiebungsvektors durch Ableitungen des Ver- 
zerrungstensors (Gl. (45)) ausgedrückt, so ergibt sich 

E",+ Te (2d,—dı |“) - Um: d+Ted)ı+tp"+P"d=0. . . . (68). 
Den so gewonnenen Beziehungen liegt die Beziehung (58) zugrunde, welche streng genommen 
heuristischen Charakter hat. Nachstehend soll zugleich eine andere Auffassung erörtert 
werden, welche ebenfalls nahe liegt. 

Betrachtet man die Komponenten des Spannungsvektors vor und nach der Deformation, 
so kann mit der Verzerrung des Materials auch eine gewisse Verzerrung des Vorspannungs- 
tensors angenommen werden. Die einfachste Vorstellung wäre jene, wonach die Komponenten 
des an der Fläche «= konst. wirkenden Spannungsvektors dieselbe Verzerrung erfahren wie 
die Komponenten des Linienelementes, d. h. 


Te ygeh:da= Ta: 2: 5 SEE. 
Diese Beziehungen verlangen wegen d#+=c dx", d#k— cd" 
Te yg= Tür yy' oder jr Ha-D) ...2.22.2. 
Dieser neue Tensor ist von vornherein symmetrisch und kann daher unmittelbar mit T*“ 
identifiziert werden. Gl. (51) liefert dann 
EUER MD 2... on nr 
und die Gleichgewichtsbedingungen (50) lauten mit Rücksicht auf Gl. (60) 


er + v";. Per. . . 2, 3 
bzw. 
Frei, Tepe +Pdi=0 ....2.2220. (6). 


Im Falle der Isotropie besitzen diese Gleichungen gewisse rechnerische Vorteile gegenüber 
den Gl. (63) und (64), wie im neunten Abschnitt gezeigt wird. Der Unterschied beider 
Gleichungssysteme erstreckt sich meist auf Glieder, die nur mit kleinen Verzerrungsanteilen in 





(56). 


ziert 


(57). 


(58) 


(59). 
y? u 


SSen- 


(60) 


(61) 


tors 


(62). 


Ver- 


(63). 


men 
tert 


ion, 
NgS- 
ten 
wie 


64). 
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Zusammenhang stehen und kaum wesentlich werden. Die Gl. (63) entsprechen formal einem 
von Biot [6] aufgestellten Gleichungssystem. Biot benutzte jedoch für den Zusammenhang 
der Ableitungen des Verzerrungstensors mit jenen des Drehungstensors nicht alle geo- 
metrischen Beziehungen, sondern nur einen Teil (jene, die aus Gl. (53) in kartesischen Ko- 
ordinaten hervorgehen, wenn zwei Indizes zusammenfallen), so daß bei seinen Endgleichungen 
noch Verzerrungs- und Drehungstensor gemeinsam auftreten. Die hier aufgestellten Gl. (63) 
und (68) stellen daher auch für kartesische Koordinaten eine neue Formulierung dar. Die 
Beziehungen (67) sind andererseits den von Trefftz [3] verwendeten Gleichungen ähnlich, 
ohne daß dort die physikalische Bedeutung in dem hier gegebenen Sinne erörtert ist. Auch 
erfolgte bei beiden Verfassern die Ableitung der Gleichgewichtsbedingungen aus dem Energie- 
prinzip, d.h. nur indirekt nach Voraussetzung des elastischen Potentials, während die Ab- 
leitung hier allein auf dem Kräftegleichgewicht fußt. Die Existenz des elastischen Potentials 
wird hier erst eine Folgerung aus dem Superpositionsprinzip und der im nächsten Abschnitt 
behandelten Proportionalität zwischen Spannungs- und Dehnungstensor. 


6. Das Elastizitätsgesetz und die Differentialgleichungen 
des elastischen Verschiebungsvektors. 

Im Rahmen der Theorie zweiter Ordnung stellt di. unmittelbar den Verzerrungstensor 
dar, welcher mit dem Spannungstensor #“ durch das lineare Elastizitätsgesetz verknüpft ist. 
Bei allgemeinem elastischen Verhalten (Anisotropie) wird diese Proportionalität durch die 
Koeffizienten E*“”? festgelegt, welche einen Tensor vierter Stufe bilden: 


a un 5 


Entsprechend der Symmetrie des Spannungs- und Dehnungstensors gilt auch für diese 
Koeffizienten die Symmetrieeigenschaft 

ara ., 2. Er re (70) 
Dem Superpositionsprinzip und damit der Existenz des elastischen Potentials entspricht 
ferner die weitere Symmetrieeigenschaft 

NE a ie rn 


Für das elastische Potential oder die Formänderungsenergie ergibt sich der Ausdruck 
g 1 .. 
A = o7 gau diu — o7 Eruve diu d, "A . . . . . . . . . (72). 
Die Gleichgewichtsbedingungen (63) bzw. (68) gehen hiermit über in folgende Differential- 
gleichungen des Verschiebungsvektors: 


Wr WR 


Hierbei ist Z?“[] ein tensorieller Operator, der bei Bezugnahme auf ül. (62) folgende Form 
hat: 

\ en 1 ER 1 u 
Tan j= (Ber + Z gen De Perg ge Pr gr Br) ler 


| 2 TE 


1 “v 2 n Au » 
+7 P+H PR- Tell 


bzw. bei Zugrundelegung von Gl. (68): 
D«[]=(Ere? + gi“ 1’) [ ve + g%” P“[],. ® i ’ . 4 , ü a (74b). 


Mit der Einführung dieses Operators ist damit das Gleichungssystem bereits formal auf das 
der anisotropen Kontinua ohne Vorspannungen zurückgeführt, und es handelt sich nunmehr 
um die Integration partieller Differentialgleichungen, deren Koeffizienten veränderlich oder 
konstant sind, je nachdem die Vorspannungen veränderlich oder stationär sind. 


7. Integration der Grundgleichungen für anisoıwropes Kontinuum 
mit stationären Vorspannungen. 
Im Falle der stationären Vorspannungen gilt 
ee Re ER 5 


Wird noch vorausgesetzt, daß anfänglich keine Massenkraft vorhanden ist und daß der durch 
die Deformation bedingte zusätzliche Massenkraftvektor p“ eine reine Trägheitskraft re- 
präsentiert, d.h, 
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t Ö E 


p - UFeTu 


(o spezifische Masse, f Zeitkoordinate), so geht Gl. (73) über in 


[P«[V.J=0. 
wobei 


MA L 0 
IBe= Beref\ee—eog"zgll 


bedeutet. Der hierbei neu eingeführte Tensor E*’"? berechnet sich zu 


Bi 1 u | Bi 
Er+r«: En E’'": — ai Tv? gg T*' ri g*: To 


bzw. £ 
RETRO an ee a 
je nachdem Gl. (74a) oder (74b) zugrunde gelegt wird‘). Hieraus folgt der Satz: Das elastische 
Verhalten eines anisotropen Kontinuums bei Vorspannungen entspricht einer Änderung der 
Anisotropie, derart daß an die Stelle des Tensors E?’"? der Tensor E*r“° tritt. 
Die Lösung der Gl. (77) stellt ein bereits von Herglotz [10] und Weierstraß [11] 
behandeltes Problem dar. In der hier gegebenen Schreibweise ergibt sich die Differential- 
gleichung der Komponenten des Verschiebungsvektors als Determinante der Operatoren 17%“ []: 


Ea 34 &iu» I]e} JJ2« Ir [V,]=0 . . . . . . . . . . (80). 
Wird eine vektorielle Spannungsfunktion X“ eingeführt, so liefert der Ansatz 
. 1 u Be 
Veoz 2 uv &a3y 11?" Ir’ [Xe] . e R . e . e ö . : (81) 
Lösungen der Gl. (77), wenn X“ der Differentialgleichung (80) genügt, welche von sechster 
Ordnung ist. Es läßt sich zeigen, daß eine der drei kartesischen Komponenten der Spannungs- 
funktion im allgemeinen zur Darstellung der Lösung ausreicht. 


8. Fortpflanzung von Störungen im anisotropen Kontinuum 
bei stationärer Vorspannung. 


Ein aufschlußreiches Beispiel für die Beeinflussung der elastischen Eigenschaften durch 
Vorspanntingon stellt das Verhalten gegenüber Störungen dar. Handelt es sich beispiels- 
weise um die Ausbreitung, einer ebenen Wellenfront, so gilt für den Verschiebungsvektor 
der Ansatz (f} ist eine beliebige Funktion des Argumentes a," + ci) 


V=f(aga+ch, Ned 


Hierbei möge a, ein Einheitsvektor sein. Die Ausbreitung der Störung ist dann durch die 
zeitliche Wanderung jener Stellen im Medium charakterisiert, für welche das Argument 
a, x* + ct den gleichen Wert besitzt; diese Stellen liegen zu gleichen Zeiten auf einer Ebene, 
deren Ausbreitungsgeschwindigkeit in Richtung der Normalen durch ce gegeben ist. Die 
Differentialgleichungen (77) liefern dann 

Ber. kod)id=0 . : ....2 22.220. (88). 
Schließlich führt die Differentialgleichung (80) auf eine Beziehung, welche eine Gleichung 
dritten Grades für c’ darstellt. Die für anisotrope Medien bekannte Erscheinung der Existenz 
von drei verschiedenen Schallgeschwindigkeiten®) wird daher bei Vorhandensein stationärer 
Vorspannungen nur insofern modifiziert, als andere Elastizitätszahlen gelten. 


9. Integration der Grundgleichungen bei stationären Vorspannungen 
für den Sonderfall der Isotropie. 


Ein Sonderfall liegt bei isotropen Medien vor. Der Tensor E’"'? hat dann die spezielle 


Form 
5) 


Eıe—@G ggg" gr+ —f'r* j P j 1 ' f x (84) 
- ) 


m 


*, Hierbei ist zu beachten, daß die Symmetrieeigenschaften (70) und (71) für den neuen Tenser Z*'""Y nicht 
mehr gelten. 
5) Sıehe z. B. Green: Cambridge Phil. Soc. Trans., Bd. 7, 1839. 
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(G Schubmodul, m Poissonsche Konstante), und der Tensor E*“'? geht mit Bezug auf 
Gl. (79b) über in 
Ber G|gt gie + gt gic+ EI ER ae 
Wie nachstehend, abgeleitet wird, läßt dieser Ausdruck eine besonders einfache Integration 
der Grundgleichungen zu. Aus den Gl. (77) und (78) ergibt sich zunächst als Differential- 
gleichung des Verschiebungsvektors: 
’ m 2 2 l o 0° 
vr ,.+ yi Fear. . 
7 m—2 + G Gar 
’ 
Es läßt sich nämlich mit Hilfe des gleichen Ansatzes, wie er vom Verfasser [12] in der ge- 
wöhnlichen dreidimensionalen Elastizitätstheorie mit Erfolg angewandt wurde, unmittelbar 
ein allgemeiner Lösungsweg finden, der zu Differentialgleichungen von nur zweiter Ordnung 


führt. Der Ansatz lautet 
Bei: 3 ee, a 


Hierin stellt F eine skalare und ®; eine vektorielle Spannungsfunktion dar. Nach Einsetzen 


in Gl. (86) ergibt sich 


2m 2 us ud Mm Alu 1 m; it u 
= F ıt® it i 44 »t% T*(— F | 


m —ä m 


Er. . 
te — Pd) = 0 | 


Es werden zur kürzeren Schreibweise zweckmäßig folgende Operatoren eingeführt: 


2 ER ern 
() ıtG E* )e— ATZA =4. 


2 m — 0 


Bar ee" I_3 


Gl. (88) geht hiermit über in 


2 7 | 1 Tyıu Alu gi 
( -„,)4® 21 —) Äh A 34 


/ 


Diese Vektorgleichung wird integrabel, wenn ®% der Differentialgleichung zweiter Ordnung 


a a: ren re 


genügt. Die Restgleichung (91) lautet nach Integration 


a1 )ir=al+e,. ar De 


m 


hierbei ist C,, ebenso wie C, Integrationskonstante. Zur weiteren Integration wird zweck- 
mäßig an Stelle von ®4 die Vektorfunktion XX eingeführt, entsprechend der Beziehung 


2 (1- .) AR sr de I Wa a 7 


Gl. (93) ist dann wiederum integrabel und führt auf die Gleichung 
F=y+X',+ C, ah. 
wobei y der homogenen Gl. (93) 
dy=V 
genügt und C, mit C, durch die Beziehung 


m —?2 


1 a‘ 
9) 9 Y > Mn ME 5. ' 
(1 „)&(6+ 1)@ T;)= 6. 
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verknüpft ist. Aus Gl. (9) und 94) folgt dann, daß X der Differentialgleichung vierter 


Ordnung 
a a a te 


genügt, sich daher sowohl aus Lösungen der Gl. (92), als auch der Gl. (96) zusammensetzt. 
Für die Bildung der Verschiebungskomponenten kommt es jedoch nur auf die ersteren an; 
denn Lösungen der Gl. (%) fallen bei Bildung von ©, entsprechend Gl. (94) wieder heraus. 
Andererseits würden sie bei F gemäß Gl. (9) nur in der skalaren Form x, auftreten, 
welche schon allgemein in der Funktion y enthalten ist. Somit kann ohne Einschränkung 
der Allgemeingültigkeit die vereinfachende Festsetzung getroffen werden, daß X; nur Lösungen 
der Gl. (92) enthält, d.h. es sei 


ee ae 5 nr re 


Diese Gleichung enthält infolge der Regeln der kovarianten Differentiation im allgemeinen 
nicht nur eine, sondern alle Komponenten X;. Daher ist es für die praktische Durchführung 
zweckmäßig, zuerst die kartesischen Komponenten X, aus der analogen Gleichung 


a er a ee DEE 2 


welche nur eine Komponente enthält, zu bestimmen und danach X, aus u=cd X; zu be- 
rechnen. Der Operator J dagegen kann im neuen System geschrieben werden. Nach Kennt- 
nis der Funktionen X, und y wird die Spannungsfunktion F gemäß Gl. (9) bestimmt, 
welche der Differentialgleichung 

fi, WER a a 


genügt. 
Die Verschiebungskomponenten errechnen sich dann aus 


/ \ 


IWw 
26 a=-Fı+2li-.)iR . ET 


m 


. 


Der Spannungstensor kann dann entweder aus den Verschiebungen an Hand der Gleichungen 


/ 9) 
Belle)... . ee 217 


welche aus (69), (42) und (84) folgen, oder direkt aus den Spannungsfunktionen, wobei die 
Beziehungen 


r 1 7 r Fu Pr v £.) ZT yı 
Re — FR +1.) Ale + 894 GP, +21) Axt.) . (104) 


m 


maßgebend sind, bestimmt werden. Die Gl. (104), welche sich aus (103) nach Einsetzen der 
Ausdrücke (102) ergeben, erlauben noch eine Kontrolle der Rechnung durch Einsetzen in 
die Gleichgewichtsbedingungen (77). Eine kurze Rechnung zeigt, daß diese in der Tat erfüllt 
sind. 
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Biegung gekrümmter, dünnwandiger Rohre*). 


Von Helmut Karl in Dresden. 


Unter Benutzung des Satzes vom Minimum der Formänderungsenergie und 
des Prinzipes von Castigliano werden obere und üntere Schranken für 
die Größe der Biegung gekrümmter dünnwandiger Rohre gewonnen. Zu- 
nächst wird das Verfahren auf schwach gekrümmte Rohre bei Biegung in 
der Krümmungsebene unter der Annahme, daß die Querkontraktionszahl 
Null ist, angewendet. Weiter wird dann der Einfluß der Querdehnung be- 
rücksichtigt, stärker gekrümmte Röhren untersucht und die Biegung quer 
zur Krümmungsebene behandelt. 


A. Die wichtigsten Bezeichnungen: 


R= Krümmungshalbmesser, 
a= Neigung eines Quer- 
schnitts gegen die lot- 
rechte Symmetrieebene, 
Halbmesser der Rohr- 
wandmittelfläche, 
-Winkel zum Punkte P 
der Rohrwandmittel- 
fläche von der Sym- 
metrieebene 2 —v ge- 
messen, 
—= Wanddicke, 
v = Querzahl, 
la,(4da)= Änderung des 
Winkels a, (da), 
i=B-tlr, 
A=R-tjr.yl—»v, 
0= r| R, 


Bild 1. 


(Bezeichnung für eine beliebige Funktion F), 


M = Lastmoment, 
M — Querbiegemoment in der Rohrwand für das Ringelement da, 
A = Formänderungsenergie, 


A = Formänderungsenergie je Raumeinheit. 


B. Biegung in der Krümmerebene. 


I. Eingrenzung der Verbiegung bei sehr großem Krümmungshalbmesser R. 


Die Wand des gekrümmten Rohres erhält aus den Biegespannungen Ablenkungskräfte, 
die eine Verformung des Kreisringquerschnitts bewirken. Hierdurch wird seine Nachgiebigkeit 
vergrößert. Die nach der Biegetheorie des geraden Rohres berechnete Verbiegung ist daher 
noch mit einem Faktor x zu multiplizieren, der eine Funktion des charakteristischen Wertes A 
ist, die wir durch Eingrenzung bestimmen wollen. 


Nach C. Weber) lassen sich die Verschiebungen von Punkten elastischer Körper, für 
die das Hookesche Gesetz gilt, durch die Anwendung der beiden Minimalsätze eingrenzen. 
Es wird dabei einmal ein Verformungszustand angenommen und die gesuchte Verschiebung 
aus dem Minimum der Formänderungsenergie berechnet (Prinzip der potentiellen Energie). 
Das andere Mal wird die Verschiebung aus einem angenommenen Spannungszustand nach 
dem Satze von Castigliano bestimmt (Prinzip von Castigliano). Ist «„ die wirkliche 
Lageänderung des betrachteten Punktes und bezeichnen «, und u. die näherungsweise be- 
rechneten Verschiebungen nach dem Prinzip der potentiellen Energie und nach dem Prinzip 
von Castigliano, so gilt 

u,.ltu SU. 

*) Auszug aus einer Dissertation an der Technischen Hochschule Dresden. Referent: Prof. Dr.-Ing. C. Weber, 
Korreferent: Prof. Dr.-Ing. K. Beyer. 

1)C. Weber: Z. angew. Math. Mech. Bd. 11 (1931), S. 244 45. 
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Der wirkliche Wert «,, ist also durch «, von unten, durch «. von oben eingegrenzt. Durch 
Verfeinerung der Ansätze verringert sich der Unterschied der beiden Näherungswerte u, 
und «. immer mehr, bis sie schließlich ineinander übergehen. 

Entsprechend den beiden Minimalprinzipien sind zwei Rechnungsgänge notwendig, um 
die Verhältniszahl x zu berechnen. Th. v. Karman?) behandelt bereits dieses Problem nach 
dem Prinzip der potentiellen Energie. Um aber die beiden verschiedenen Verfahren besser 


miteinander vergleichen zu können, wird hier die Karmansche Untersuchung mit einem all- 


gemeinen Reihenansatz für die Querverformung gebracht. 

Das Prinzip von Castigliano liefert eine obere Schranke. Wir nehmen am Rohr- 
querschnitt Längsspannungen an, die mit den äußeren Kräften im Gleichgewicht stehen. 
Hierbei tritt in der Rohrwandung eine Querbiegung ein.: Der Faktor x, der hier aus der 
Spannungsenergie bestimmt wird, ist die obere Eingrenzung der wirklichen x-Funktion ’°). 

Wir beschränken uns zunächst auf schwach gekrümmte Rohre und wenden das Ritz- 
sche Näherungsverfahren an. Bei der Eingrenzung mit Hilfe der beiden Minimalsätze wird 
der Halbmesser r wie auch die Ordinate y gegen R als klein vernachlässigt. Das Krüm- 
mungsverhältnis o wird mit genügender Genauigkeit gleich Null gesetzt. Die Querdehnung 
bleibt zunächst unberücksichtigt, die Zahl » wird gleich Null angenommen. Der Einfluß von v 
und o auf die «-Funktion wird erst später untersucht. 


a) Prinzip der potentiellen Energie. 


Formänderungsansätze, Formänderungsenergie. Wir betrachten ein Rohr- 
teilchen von differentialer Länge, das durch zwei Schnittebenen abgegrenzt wird, welche durch 
die Krümmungsachse k—h gehen und unter dem Winkel da gegeneinander geneigt sind. 
Wir erhalten ein Ringelement. An diesem werden nach Bild 2 folgende Formänderungen 
angenommen: 


Bild 2 
(links). 





Bild 3. 


. Die Änderung Ada des Zentriwinkels da durch Verdrehung der Endquerschnitte. 

2. Die Querverformung des Kreisringquerschnittes in ein Oval. Nach Bild 3 ergibt 
sich eine Verschiebung w, die in die Komponenten mw, und ®; zerlegt wird. (w, 
werde nach außen, w; in Richtung wachsendes » positiv gerechnet.) Die Vertikal- 
komponente von w ist my. 


Die Längsdehnung e, besteht aus zwei Anteilen, die von der Änderung des Zentriwinkels 
und der Querverformung herrühren. 


—_#_4da_ r Ade RE... SA. 
AT y+-Rda Rda "Pp ig 
g=&-+Ee. 
Hierbei wird vorausgesetzt, daß die Rohrquerschnitte eben bleiben. Im Nenner wird y gegen 
R als klein vernachlässigt. 
Für die Querbiegung infolge der Querverformung ergibt sich aus den Verschiebungs- 
komponenten w,, m; die senkrechte Durchbiegung 





wy=W,C0OSP — nı-SINQ. 
?2) Th. v. Karman: Z. VDI Bd. 55-(1911), S. 1889. ’ 
»),C. Weber hat sowohl für die Biegung in der Krümınerebene, als auch für die Biegung senkrecht hierzu, 
die später behandelt wird. die Grundlösung nach dem Prinzip von Castigliano als obere Grenze gefunden. Hier- 
über besteht jedoch keine Veröffentlichung. 
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I/m eine Beziehung zwischen den beiden Funktionen w, und w; zu erhalten, wird die Längs- 
dehnung &„ im Ringelement gleich Null gesetzt. 


1 
EmMH= —_ ww, + w;) =», 
’ (Ableitungen nach x sind 


w,=——n,, durch ’ bezeichnet). 
wy„= — (m - sing + m;-cosgp) 
Hiermit: 
1 Ida 2 - M 
eg, =7!|f- cos ( n'; SIn < - mW, COS ( ne er Fe s 
1 R da / t 7 t Y 


Für das Biegemoment M im Querschnitt des Kreisringes gilt: 


M 1 M 
ir. (w, + w,) | 


BJ 

oder (2). 
wi | 
Ejmre) 


Unbekannt ist die Funktion w;. Sie muß so bestimmt werden, daß die Formänderungsenergie 
zum Minimum wird. Die Formänderungsenergie wird hierbei durch e, und M, also durch w; 
ausgedrückt. An Stelle von »; führen wir die Funktion 


IR da si EN ) 

= 74, Mr sing ee a a 
ein. In G]. (1) und (2) eingesetzt erhalten wir 

ER 

AR da 

M Ada ( WY 

EJ rda En 


5 Een 
q 
(5). 


Nun verwenden wir für ®»; den allgemeinen sin - Ansatz 


n—=.—x 
r-dda \” 
4=2- / sin®2ny. 
n,: —=2 2 /, "sin2ny 
> 
Hieraus folgt: 
n==x 
hj - 
B—- ) [En +VD)cos@n -Yp+l2n—1)cos(n+1)g], 
Ey 
N N=—= X 
ap v en , 
( je- Y 4n(4n?—1)-c„-cos2np. 
cosq du 
ee} 


Die: zunächst unbekannten Faktoren c„ ergeben sich aus der Bedingung, daß die Form- 
änderungsenergie zum Minimum werden muß. Die Formänderungsenergie wird: 


Erce Im 
A=5 \\le.ar+sin,-as 


7 - 


—_ Er t (Ada) 2 | [eos B’dp+z | 





EJEGE 


"Em “8 
0 0 
_Er’t (Ada)? ( - er res ae er u 
un: De 7° 12 | je0s7 I enlten+1)cos(n )p+l@n-1)eos(@n+1)p]| dgq 
) 1 
(6) 
82 f > 2 
+ \ | > C„n(4n?—1)cos2n 2 del 
u l 
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Bedingungsgleichungen. Aus dem Satze vom Minimum der Formänderungs- 


energie folgt: 


A 3 z 
= -—2\| (cos > [(2r+1l)es@n-—1)p+(2n- cos @n+ gl) 
v 1 


(2u +lD)eos(2u -1N)p+(l2u—1N)eos(2u+1l)g]ldg 
sat 
3 | 
v 1 
Rechnen wir diesen Ausdruck aus, so erhalten wir Gleichungen folgender Art: 


. 


en n(4n?—1)cos2np-u-du’—1)cos2ugp-dp. 


Aus, tur: C,+ EERRY +aug CR =bu. 


Für den vollständigen Reihenansatz ist k= zu nehmen. Durch Abbrechen nach dem ersten 


oder einem folgenden Gliede (k=1, k=2 usw.) ergeben sich Näherungslösungen, deren Ge- 


nauigkeit mit der Zahl der Unbekannten wächst. Für die Vorzahlen a.„ und die unabhängigen 
Glieder b. gilt: 
20, =2\[@n +1) ceos®@n —N)p+l2n—1)cos(2n+1)y] 
0 
(24a +l)cos2u -Dp+l@u—1)cosQu+l)oldyg 
gj: > 
+ z n» u (4n’— 1)(4u® —1) \ cos2ngy-cos2ugp-dg, 


n-bu=2 \ cosp[Z2u+1)cos(2u— 1 +l2u-—-VDeos@u+Vg]ldp. 
0 » 
Hierin sind « und n ganze Zahlen. Durch Ausintegrieren erhält man 
4n? 
Ayun=2(4n’+1)+ 3 (An? - 1)?. 42 

An-nn=(2n—3)(2n+1) Ann en —1)(2n +3) 
Alle übrigen a.„ =0 
=, 
Alle übrigen b=0. 


Für die unabhängigen Werte und Vorzahlen erhalten wir nachstehende Matrix: 





1 2 3 4 b 
1 2(5-+62°) 5 3 
. 1 
2 5 2 (17 + 600 4°) 21 
3 21 2 (37 + 7350 4?) 45 
4 45 2(65+ 423367? 














Berechnung der Faktoren x. Der unter (6) angegebene Ausdruck für die Form- 
änderungsenergie kann für die Bestimmung der x-Funktion durch die Bedinguugsgleichungen 
vereinfacht werden. Zunächst wählen wir eine vereinfachte Schreibweise für Gl. (6) 


7 


9312 [ 
(G- Eon H’dp+ 3 | len In? dg 


v 


Er’t (Ada) l 


ou 2R da 


Sa, 


G, H„, J„ sind Funktionen von g. 
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Die Minimalbedingung lautet allgemein: 
0A 7 


; ö BE : . 
3, = —2|(@ - 206. H„)H,-dg ” 3 \r Leu In dp. 


0 0 
Mit c.==0 kann sie auch geschrieben werden: 


7 r 


24° b 
0 2c; \ (G 2 c„. A, Hı«-dp+ 3 cH\ IR: Lem Judy . 


0 v 
Die Summe aller dieser Ausdrücke ist 
Ä Er IT 
0=—2|\ Fe). H,(G— Ic„H,)dp+ 3 te WE & 7 77 Fu 
0 0 

Nach Subtraktion der Gl. (8) von (7) erhalten wir 
Er’t (Ada) i 
— . 2G(G — Zc,: Be BEN 9 
IR 5 \26 C„:H,)de . \ (9). 


A: 


Daraus folgt mit 


G=c0sy, H„=(2n + 1)cos(2n -Dpo+(2n l)cos(2n+1)y, 


Er’t; .M 
A= ng. (Ada’)-[1-30)]: >= Ida, 
fd _MRda 1 
Ca> Ertai -3c, ' 
nn - (10) 
1—3c, 


Wir berechnen nun die x- Werte, indem wir den Reihenfaktor ec, bestimmen und in (10) 
einsetzen. Für c, bestehen verschiedene Näherungen, je nach der Zahl der Gleichungen und 
der Unbekannten, die zur Bestimmung verwendet werden. Obenstehende Matrix zeigt, daß 
die Gleichungen aus dreigliedrigen Ausdrücken mit Symmetrie zur Hauptdiagonale bestehen. 
Bei Entwicklung von x in einen Kettenbruch*, wird 


1 
er 
35, 
1 
a. 
d,; = 
a ad, 
1. Näuerung ai s 
N ee AU Bi 


33 


2. Näherung 


3. Näherung 


Hieraus: 


a A _10+12# 
1. Näherung: x,= irB® 


105 + 4136 4? + 4800 4! 


> 6° > . e- ws . , = Ar 
2. Näherung: %,= 3 556 774800% 


2 äes 2... 82473912 4? + 24461764* + 2822400 4° 
» NAherung: #= "3, 32804? + 3293764 + 28224004 


Statik im Eisenbetonbau, ?2. Aufl., Bd. I, Berlin 1933, S. 230. 
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b) Prinzip von Castigliano. 


Spannungsansatz und Querbiegung. Für 
die Längsspannungen wählen wir den Ansatz: 





M Fe 
o0=—, cosyp-+ \ ee. (11). 
r” tn| — l 


1 


Sie erzeugen in den Längsfasern des Ringelements mit 
dem Zentriwinkel da Ablenkungskräfte o-dF-da 
(Bild 4), die sich über den Umfang des Ringelements als 
Streckenlast verteilen. Diese ist durch folgende Gleichung 
bestimmt: 

a=o-.t-da. 


r 
Die Ablenkungsbelastung ruft im dreifach statisch unbe- 
stimmten Ringelement ein Biegemoment hervor. Für 








dieses erhalten wir Bild 4. 
Mda j = Cn | 
S — — inc 9 B. " | os d .ene > i ‘ 
M 17 [cs 24 2 aa +D@n rm ®tDVecos2ng+n cos? (n+D) ol; Ai 
1 


Spannungsenergie unter Berücksichtigung der Minimalbedingung. Die 
Spannungsenergie wird: 
1 Bu a a 
A=55\\ \"dv+\anyds. 
Setzen wir für o und WM Gl. (11) und (12) ein, erhalten wir in vereinfachter Schreibweise 


be) T & 








M’Rda R u I’ 3 f ET > 
Er’tm? \ G+ 5. eu Hu] de+tzz| J+ I On Kl d P\ 40. 
v 1 0 1 


Hierin sind @, H„, J, K„ Funktionen von 9. 


Die Minimälbedingung ergibt: 


0A 3 er , 

zu. 6+ Ei RR Hu: H,.:dg +7 „+2 : eu: Ku Ky:do a 
0 l 

Setzen wir die Summe aller ae in den ae (13) ein, so erhalten 


wir ähnlich wie beim Minimalansatz der potentiellen BER 











M’Rda : 2 B 
 &_ — E r?’ t n? \ G G Zu Bl Ca ® 1, d q +4 J u >& ° u Ku) d q 
Mg 1 
und nach Integration: 
. M Rd af 3 je =, 4 15 
£ "»BErtn| +752 +7 > AI E » ° ö r . . R ; (10), 
Hieraus finden wir 
(d Hei c,\ 
das gytn erg 
ur 3 ce, / * 
»<1I+,5t7 j? P . . ; . . s . A . ° s ; " . . (16). 





Die «-Funktion enthält den Teil 1+3/4-4°, den wir als Grundlösung bezeichnen, und den 
vom Näherungsansatz abhängigen Teil e,/4 4°. 
Berechnung der Faktoren x. Aus den Bedingungsgleichungen (14) erhalten wir 
ein Gleichungssystem von der Form 
AnıC, + Aue, + ER + Au: c; = bu. 


Durch Integration finden wir: 





ö 3[n?+(n + 1)?] 
an U ER R 
ER rom (n+)(2r +1] 
3 3 
n-nn =Iyrldn? 1)’ Haan tiPR@n+Dent+E 





a 
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Alle übrigen aun =. 1 
I, .) 


Alle übrigen b.=0. 
Diese Werte sind, multipliziert mit 600, für n=1 bis n=4 in der nachstehenden Matrix 
zusammengefaßt. 

















1 2 3 \ b 
1 25 (487? +5) 15 300 
. : 20 
2 15 1200 4? + 13 Be 
‘ 
20 25 j 25 
3 = : (9408 4? + 25 
7 196 28 
>57 
_ 13200 4?+ 41) 
28 36 
Nach den früheren Ausführungen läßt sich ec, als Kettenbruch darstellen. Setzt man ce, 





2). in (16) ein, so erhält man: 
3 
Grundlösung: #,=1 , 
)ie , re ” 44° 
* 3 3 
1. Näherung: „14757 page HB) 
3 3600 2? +39 
> 5 3 + a Pe a . z — 
” 2. Näherung: »,=14 74: 57600 + 66410 +56 
2 8 11289600 ä* + 152304 42 + 261 
3. Näherung: 1, = 14 775 7 180633600 4° £ 21252864 2° + 320792 7 + 360 2° 
3). 


II. Berücksichtigung der Querdehnung. 

Die Ableitung der Ergebnisse der vorhergehenden Abschnitte wurde durch die Annahme 
vereinfacht, daß die Querzahl »=0 ist. Ergänzend wollen wir nun den Einfluß der (Quer- 
dehnung auf die x- Funktion feststellen. Hierbei benutzen wir das Prinzip von Castigliano. 

Für unsere 
}). Untersuchungen 
können wir wie 
früher annehmen, 
daß die Spannun- 
gen senkrecht zur 
Rohrwandung o, 
(Bild 5) keine Be- 
deutung besitzen 
und setzen sie 
‘ deshali» gleich 
Null. Fs ıst dann 
ein zweiachsiger 
Spannungszustand 
vorhanden. 

Die Längs- 
spannungen des 
Rohrkrümmers 

2) nennen wir jetzt .. 
%,, während wir die Spannungen aus dem Quer- 
biegemoment M mit 0%, bezeichnen. Letztere 
haben infolge der Querdehnung das Bestreben, 
den Kreisringquerschnitt zu verwölben. Da dieser 
aber eben bleiben muß, müssen in Längsrichtung 
neben den früheren Biegespannungen o noch 
Spannungen wirken, welche die Verzerrung auf- 
heben. Wir setzen also: 






N 
| 
| 
I 
I 
| 


Näherung potentielle Energie 





ZZ 


Castighiano 


”„ ” 


rundlösung Castigliano 


un AO NR 
SWS Sum» 





| 
| 
I 
| 
L-Aose] Bild 5. 











Asymptote für A- 


nn m 1 m {17 -  —  _ —- 








0 2 098 9 05 06 07 08 039 7 
Z MON A 
Er. . 50%. (17). Bild 6. Eingrenzung des Faktors «. 
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Wir verwenden für o den früheren Ansatz (11). Die zusätzlichen Spannungen vo, rufen 
Ablenkungskräfte hervor, die sich jedoch in jedem Querschnitt irda des Ringelements 
selbst ausgleichen, so daß die Ablenkungsbelastung gegenüber dem früheren Ergebnis keine 
Änderung erfährt. Da die Spannung v-0o%, keinen Einfluß auf die Querbiegespannungen aus- 
übt, erhalten wir für das Querbiegemoment W den früheren Ausdruck (12) und hieraus o,. 


Die spezifische Spannungsenergie für den zweiachsigen Spannungszustand wird 
er 1 
A=;5 Rt — 270. a re Yen er t a W 
Setzen wir 0.=0+1v:0, ein, so erhalten wir: 
A= Lord — v?) 07] (19 
=5% o EP a ee EL er J), 


Bei Integration von A wird: 


y—=}!ın 
pe 
-dV+A-r) \opyde: 
Pe g=0 


Für o und M sind die früheren Ansätze (11) und (12) einzusetzen. Hieraus folgt, daß bei Be- 
rücksichtigung der Querdehnung nur der charakteristische Wert A in den «-Funktionen nach 


an j 4 
dem Prinzip von Castigliano in A= 77 übergeht, während die funktionalen Be- 


ziehungen ungeändert bleiben. Durch fortgesetzte Steigerung der Genauigkeit der Näherungs- 
ansätze erzielen wir auch hier Ergebnisse, die in den wirklichen Wert von x übergehen. Eine 
Untersuchung nach dem Prinzip der potentiellen Energie erübrigt sich daher. 


Setzt man sowohl in den drei Näherungen nach dem Satze der potentiellen Energie als 
auch in der Grundlösung und den drei Näherungen nach dem Satze von Castigliano A an 
Stelle von A, berechnet dann für verschiedene A-Werte den Faktor x und trägt ihn über A 
auf, so erhält man Bild 6. 


In dieser Darstellung läßt sich die schrittweise Steigerung der Genauigkeit der Ein- 
grenzung gut verfolgen. Für größere A-Werte ist die Grundlösung des Prinzips von Casti- 
gliano eine genügend genaue Näherung. Legt man die für technische Berechnungen meist 
ausreichende Rechenschiebergenauigkeit zugrunde, gilt folgendes: Nach dem Satze der poten- 
tiellen Energie stimmt die erste Näherung bis A=0,6, die zweite bis A=0,2 und die dritte 
bis 0,1 > A >0,05 mit den höheren Näherungen überein. Beim Prinzip von Castigliano ist 
dies für die erste Näherung bis A=0,35, für die zweite bis A=0,10 der Fall. Der Genauig- 
keitsbereich der dritten Näherung läßt sich nicht festlegen, da die Vergleichsmöglichkeiten 
fehlen. Aus der nachstehenden Zusammenstellung geht jedoch hervor, daß für A=0,05 der 
Fehler bei Verwendung der dritten Näherungen beider Minimalansätze kleiner als 2,2 °/o sein 
muß. Also sind hierfür beide noch brauchbar. Kleinere A-Werte als 0,05 werden praktisch 
kaum vorkommen. Für die dritten Näherungen erhält man: 


A= 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 


Pot. E.: 33,991 17,307 11,555 8,6751 6,9223 5,7280 4,8598 4,2086, 
Cast.: 34,752 17,316 11,555 8,6751 6,9223 5,7280 4,8598 4,2037. 


A=-05 05 06 07 08 09 Te 
Pot. E.: 3,6948 3,2929 2,70% 2,3158 2,0406 1,8415 1,6984, 
Cast.: 3,6949 3,2929 2,70% 2,3159 2,0406 1,8415 1,6934. 


III. Berücksichtigung der Krümmung der Rohrmittellinie. 


In den bisherigen Untersuchungen wurde angenommen, daß der Krümmungshalbmesser R 
gegenüber dem Halbmesser r des Rohrquerschnitts sehr groß ist. Es soll nun festgestellt 
werden, welchen Einfluß die Krümmung des Rohres auf die vorher behandelte Formänderung 
hat. Zu dem bisherigen «- Wert tritt noch eine Funktion von o und A. Diese gibt den An- 
teil der Krümmung an dem Faktor x an, der das Verhältnis der Verbiegung des gekrümmten 
Rohres zum geraden Rohre darstellt. Sie wird nur nach dem Satze von Castigliano be- 
stimmt, da eine genauere Annahme über den Verformungszustand zu kompliziert wird. In 
den folgenden Untersuchungen werden nur die Glieder mit o und o? berücksichtigt. 











En 


|: 


ans 
zus 
Fal 


ode 


Pot 


nüg 
der 
Ver 
wir 
Au: 


Die 


gli 


M, 








). 


Lo Zn) > BE ZW 











Z. angew. Math. Mech. 


Bd. 23 Nr.6 Dez. 1943 Karl, Biegung gekrümmter, dünnwandiger Rohre 339 





Der Gang der Berechnung ist der gleiche wie früher, nur muß jetzt der Spannungs- 
ansatz geändert werden. Wir bilden ihn aus dem bekannten Ausdruck für den Spannungs- 
zustand des gekrümmten Trägers, den wir durch eine cos-Reihe mit zunächst noch unbekannten 
Faktoren verbessern. Die Freiwerte erhalten wir aus dem Minimum der Spannungsenergie. 


Spannungsansatz, Querbiegung. Längsspannung: 


I. 
M [lo 4 — 30°? cos g \ v 
= 5+ .— — c„-cos(n+1)g j (20) 
rin! t 1-+ocosg 7ER 


1 
1 
1l+ocosgq 


oder nach Entwickeln von in eine Reihe unter Vernachlässigung aller höheren 


Potenzen als o*: 


[1 
fm | P— 5eos2y+ 4 cos3r + „cos(r —+1)g i . (21) 
r. IT | . 


l 


x 
M | o ” F 
— 
l 


Der schwach trapezförmige Querschnitt der. Wandung des Ringelements wird mit ge- 
nügender Genauigkeit als Rechteck angenommen. Dann ist für parabelförmige Verteilung 
der Schubspannung die Schubverteilungszahl 1,2 zu nehmen. Die Querzahl wird wieder zur 
Vereinfachung der Berechnung gleich Null gesetzt. Unter diesen Voraussetzungen erhalten 
wir für die Schnittkräfte aus den Ablenkungskräften im statisch unbestimmten Ring folgende 
Ausdrücke: 


Biegemoment: 





Mdajo 12 +0? 0 IR; 
M: im |9 608 P + 75 cos2p—  cosdg+ 5760844 
c | 2#+1 Bu - 
+ > [o®- + cos g +c0s3p)+ 2 (06089 +2cos?2p-+costy)} ;° 
3\0 79 ro“ 
+2 | > BORN +5c0s3p +3cosdg & \ 
"30 48 77 SE 
Längskraft: 
Mdaj| c 
N = 5-— 11—cos2 > (2cos cosdY 
ak 7+ 5 (2cosg 0s3p) 
(23). 
C, > N cz D - 
77 (cosög -cos4p)+7 (cosdg c0S5Y)..... N 
Querkraft: 
= Mdaf. Be -; c - a 
V=5 .= Isin 29 +5 sinsp+ 3 (sin 2p7+sindy) 
(24). 


6 ee m 
+7 (sin3g- sindp)..... 


Spannungsenergie unter Berücksichtigung der Minimalbedingungen. 
Die unbekannten Faktoren c,„ des Spannungsansatzes werden nach dem Satze von Casti- 
gliano ermittelt. 


p=?27 y=’# y=2a 9y=2ı 
7 1 ( M cp 12 [ ©° 
A=yy \ ®dv+ >E \ jy dst35 \ pAst+ 5 \ 48: 
r : 0 gi ns ; x 0 me 


Ein: { Ic bedeutet die innerhalb dieser Klammern stehenden Ausdrücke von 





MM’ Rdajtl > BL 
A = Eritr \ cosp+ I en cos(n Hl) p| (I+ocosp)dg 
0 l 
ie de ei el) de \ 
> E% er ge 4. > . 
ja \ \Iml+ocosy 7 \ fe 1+0c0osy ”.. \ \joaıl + 0ocosg]| h 
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Wir wählen wie früher die vereinfachte Schreibweise mit @, H„, J, K„, L, On, P, Q. als 
Funktionen von g. 
M’Rda 1 | g 


A= a Tat \IE+ Fon HP + eospdp+3\[J+ Le. Au] 


7 7 


dgq 


h 1+0ocos Q 


0 0 
7 7 
j . BR: i 
+e#\IE+ Fe 0]. 1+0cosg +0°: 5 \IP+ Zen Qu] 


( d Yo d qQ | 


1+ocosgf' 


T 7 


=0=4#? | Hx-[@+ ce, H„](1+0c0sp)dp+3 | Ku [J+Fc,K 


v v0 


0A dy 


OC; 


n] 1l+0c0sg 


2 j2 N dp jeir, ,„12# ! ah, ; 
#4 \ 0,.[L+2e, O,] I+oc0s9 0 3 \ Qx;[F +2Cc, Qu] f 
v v 


{ dgq 
+ocosp' 


Durch Ausrechnen erhalten wir eine lineare Gleichung folgender Form: 


bu=Mu:C5, tUu:C5 +. Auk'Ck. 


BERN: I da 
Hierin ist nach Division von = durch x: 


7 7 


44 3 
aun= | H„Hull+ocosg)dp+ \ Kuku 


dg 
1 = 0 COos q 


0 1) 


DL dgq 1242 [ dgq 
S ( ( v l P. S m 
+ \ LT erwreee der \ 99 a 
v M 
j44? E.: =: a dg 
um — GH. °0S 9 /) J-Ku° 
b | 2 \ H (l+ocosp)dp+z \ J-R IF ocosg 
1) 0 
2 7? e d op b 12 A? ’ d op | 
- 4 RR Pr = < "- . 
Baer" \ r I+ocosp"* 5n \} bi l+ocosgl 


u 0 
Unter Berücksichtung der Tatsache, daß c, und ce, die ungeradzahligen Potenzen o, o° 
enthalten, ergibt sich für die Bedingungsgleichungen das Schema Seite 341. 
Wird die Berechnung der Spannungsenergie wie früher mit Hilfe der Bedingungs- 
gleichungen vereinfacht, ergibt sich: 


7 7 





M’Rdafl 5 -{ dgq 
sar!e. 2. ENTE PR . 
“ Erin \ sw +2. Aujtitee At dr2e u. Er 
0 u 
vo? ( dg se( dgq | 
+7 4 4 & ß + ® ’—+ 2 , 
iR \ ib Lt OL 90089 > \ ee en u) oc0sg] 
0 u 
oder ausführlich geschrieben 
M’Rdajl 0 0? 2 o,+-o? 
we 777 \\ cosy — > cos2gq +7 60534 .($+ { cos p)dg 
0 


0° c ©; 
+ T \ ji cos2p)+ > (2cosgq — cos3g)+ 3 (cos2gp —cos4y).. 1 (1 —cos2y)dy 





i REN RE Be : 
sin?2y+zsin3y +7 (sin2p +sin3g)...}sin2y -dy 
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Grundlösung 


Näherung 


Näherung 


Grundlösung 


Näherune 


. 


2. Näherung 








> 0? ’ . 
120 + e (1134? + 145) 





40 (124? + 1) 


2 
30 0 (842 — 1) | 10(484? +5) + 5 (272 2 + 75) 


60 (40 4? — 3) 


60 (404? 


: (600 A? 


60 (404? 


2 


60 (4032 


x 0? a 6 
(1200 4? + 13) + —- (8167? + 95) 


25 7 


16 


0 (2864% — 3) 45 
249 


o (280 A? 


(7350 A? + 
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- 


3 (fo 12-+0°? 0 i ME: BU ; 

+5 \1$ cosp+ CH cos2Y - g 6053 +54 c0s4gq + 30’ 3 608g +cos3p) 
0 

C, c; = . 5 
Fr 6 (ocosp +2cos2p + cos4tyg )+3 7 9°cosg +5cos3p + 3cosäg 

> «€ 

C 12 +50? 20 30° 
= 3,8 cos4o +2cos6p).. N [ E 00829 — 5 c0s3p + g cos 4 pP) dg | i 


Durch Integration finden wir mit 
oh, s=Kk+ef:k, . 
D 





,=ok, co=k,+o’k.. 
M’Rdaj| 3 1 
ii“ a 
=Sprialtrgetrgahte | 
. + aM. 
132045 1 Be u 
ee a ze ha + aa all 


Berechnung der Faktoren x. Durch Ausrechnen der Reihenfaktoren mit Hilfe der 
Matrix erhalten wir für das Prinzip von Castigliano je nach der Zahl der Unbekannten, 
die wir berücksichtigen, verschiedene Näherungen. des Korrekturgliedes, die zu einem ent- 
sprechenden #- Wert für sehr großen Krümmungshalbmesser R hinzutreten. 


Grundlösung: 
2 11 1 
me (tan gear) mh, 
1. Näherung: 
11 65120 A +2U6I6TS AA +SOMSA HS AHBLTE) )_y 





en = [11 4 344: = 11090900 A + BISHER + BSEIDT ICH 148002" TS 
Die A-Funktion des Korrekturgliedes der 2. Näherung ist bereits ein sehr umfangreicher Aus- 
druck. Daher wird die #xx,-Funktion (n=5) punktweise für einzelne ausgezeichnete /-Werte 
bestimmt. 

Aus Bild 7 ist die Konvergenz der Näherungen des Korrekturgliedes zu erkennen. Bild 8 
zeigt die Änderung der x-Funktion bei Berücksichtigung der Krümmung. Hier ergibt sich 
ein nur unbedeutender Fehler, wenn wir in der unter der Annahme »—=0 berechneten Ver- 
besserung an Stelle von A die Größe A=Ä/yl—v” setzen. Die #-Werte sind für o=0,W%, 
0,25, 0,5 und 1,0 aufgetragen. Das Schaubild zeigt deutlich, daß die Korrekturen füro=1/3, 
die wir mit dem verwendeten Berechnungsverfahren genau bestimmen können, das Ergebnis 
für o=0 nur unwesentlich ändern. Die verbesserten x- Funktionen sind ebenfalls für {/r= 0,15, 
0,20, 0,25 und 0,30 in Bild 8 eingetragen, um den Einfluß der Wandstärke zu zeigen. 


C. Biegung quer zur Krümmerebene. 


In den vorhergehenden Abschnitten ist die Formänderung eines dünnwandigen, ge- 
krümmten Rohres bei Biegung in der Krümmungsebene untersucht worden. Nun soll die 
Formänderung durch Momente ermittelt werden, die senkrecht zu dieser Ebene wirken. Um 
von deren Spannungs- und Formänderungszustand ein klares Bild zu bekommen, wird das 
Verhalten des Rohres durch ähnliche Methoden festgestellt, wie sie in den vorstehenden Ab- 
schnitten angewendet wurden. Hierbei wird das Verfahren von Castigliano benutzt, 
da es hier einfacher ist, den Beanspruchungs- und Formänderungszustand aus einem Span- 
nungsansatz aufzubauen, als den Formänderungszustand auf dem Wege der Anschauung zu 
finden. Für die folgenden Untersuchungen wird der Belastungfall Bild 9 gewählt. 


I. Kräfte am Ringelement. 

Wie früher betrachten wir ein Rohrteilchen von differentialer Länge, das durch zwei 
Schnittebenen abgegrenzt wird, welche durch die Krümmungsachse k—h gehen und unter 
dem Winkel da gegeneinander geneigt sind. Jeder Querschnitt wird auf Biegung und Torsion 
beansprucht. Es ergeben sich die Komponenten, die in Bild 11 aufgezeichnet sind. Dieser 
Belastungszustand läßt sich nach Bild 12 in zwei Gleichgewichtsgruppen I und II aufspalten, 
wenn man die mit sina und cosa multiplizierten Teile trennt. In der Belastungsgruppe I 
tritt zu den Biegemomenten noch der differentiale Betrag aus der Änderung der Torsion, um 
das Gleichgewicht zu erhalten. Ebenso tritt in der Belastungsgruppe II zu dem Torsions- 
moment der differentiale Betrag aus der Änderung des Biegemoments. 
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Bild 8. Krümmungskorrektur x nach Castigliano 
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II. Beziehungen zwıschen Torsions- und Biegespannungen. 


Die Belastungsgruppe II ergibt neben den Torsionsspannungen eine Belastung des Ring- 
elements in Richtung der Rohrwandung. Da die Querschnitte mit der Höhe (R+r-cosp)da 
keine Biegemomente aufnehmen können, erhalten wir eine Bedingung, aus der wir eine Be- 
ziehung zwischen Torsions- und Biegespannungen ableiten können. In Bild 13 sind die Kräfte 
angegeben, die am Volumenelement mit der Dicke t bei Belastung durch Gruppe 11 auftreten. 
Wir bestimmen zunächst die Resultierende in der a-Richtung, wobei wir dr=tr’-dy setzen. 














j do 
P,„=-?trrtdpdasinpg+r(R+rcosp)tdpda-— en rtdoda. 
C 
Die Längsbelastung der Rohrwandung als £ p 
Streckenlast erhalten wir durch Division T._ /N 
. 
von P, mit r-dg ; 
Rida Pi 
Pa = Es 
F 6 x Ne 
(26). / Pie \ 
we ! e do / rw 4 
T +olt c0Sp — 2Tsin@ — I— Fest # 
| da, o_ 
I , 
Um die Beziehungen zu finden, die zwischen | Pr 
oundr herrschen, betrachten wir die Schnitt- 7 
kräfte an dem Volumenelement aus der Be- #— - ee 


lastung p.. Als Gleichgewichtsbedingungen 
erhalten wir nach Bild 14 








diM=M-de .....M, 
dM,=—-Ma:de+0Qrdyp . (88), 
Bene. = 

Da das Biegemoment M,=0 ist, folgt aus Bild 14. 


(27) M,=0 und damit M,=const. Da- 
nach wird aber nach (28) Q=const. Dies bedeutet jedoch, daß nach (29) p„ = ist. 
Die Bedingung p„=0 läßt sich mit Hilfe von (26) folgendermaßen anschreiben: 


2 , ' 00 c 
tolt cos op - ZTsinp—,, Ben et 


Bei schwach gekrümmten Rohren wird wiederum o=0 gesetzt. In diesem Falle gilt dann: 
ee A a a ENT, En: 


also r unabhängig von @. Die Schubspannungen sind damit konstant über den Kreisring- 
querschnitt verteilt und unabhängig von den Längsspannungen aus der Biegung. Berück- 
sichtigen wir o, so muß (30) voll erfüllt werden. 


Ill. Berechnung der Endverdrehung für e=0 nach dem 
Prinzip von Castigliano. 


Wir wollen die Untersuchung für die Annahme durchführen, daß die Rohrachse nur 
schwach gekrümmt ist und o mit genügender Genauigkeit gleich Null gesetzt werden kann. 
1. Spannungsansätze. 

Es werden folgende Ansätze gewählt: 
Für die Torsionsspannungen : 





an 
Msin.a = : 
Ten lt R, t„.cosn e) Be ee RE EPG ES 
1 
Für die Biegespannungen: 
Mcosa| . — u Pr 
=, |nP E= Y Cy-sın mp (33). 








1 





Die Gleichgewichtsbedingung | ) ırdF=M-sina und \ | or-sing-dF=M-cosa sind durch 
die Ansätze erfüllt. Für o=0 folgt aus (31), daß alle t„ gleich Null sind. Die Biegespannung 
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enthält in diesem Falle nur ungerade m, so daß 'wir, ohne einen Fehler zu begehen, die An- 
sätze folgendermaßen abändern können: 


Msına 
T . (34), 
2r’in 
Mecosa,. 5 : Pe 
0: 23 [sing +Z2e,-sin@n+VD)o]). . -» » 2 .2.....(8). 
r’tn 


Die Freiwerte werden wieder nach dem Satze von Castigliano bestimmt. 
Die Spannungsenergie setzt sich aus zwei Anteilen zusammen: 
a) der Spannungsenergie aus der Biegebeanspruchung (mit Querbiegung), 
b) der Spannungsenergie aus der Torsion. 
Wir ermitteln zunächst die Querbiegung des Ringelements, um dann die Formänderungsarbeit 
der Lastgruppe I zu berechnen. Dieser Teil genügt bereits zur Bestimmung der Faktoren c,, 
da diese Freiwerte im Ansatz der Torsionsspannungen nicht enthalten sind. 






2. Querbiegung. 

Belastung: An jedem differentialen Teilchen des Ring- 7 dr 

. « 3 .. 35 da 

elements greifen die Lastkomponenten P, aus der Änderung des * 9% 
Torsionsmomentes und P,„ aus der Abtriebsbelastung der Längs- 
spannungen o an. Wir erhalten aus diesen äußeren Kräften nach 
Division mit rd folgende über den Ringumfang verteilte 
Streckenlasten 


4 
Mecosa-t-da_ Mcosa 2 
Py = 2 r2.t ar 7 y 
2r?.t-n er! n Ra _L 
= De ne 
M cosa ' iu i 
pr =0o.t.da= 2 da|sinp+ > c„sin@n+1l)p|. 
’ — / 
Die Belastung ist antimetrisch zur Achse » — v. 
Querbiegemoment: Für das Biegemoment des statisch U ARE 
unbestimmten Ringes erhalten wir: na 
Bild 15. 
Mcosadaf. ce“ Ca h : | 
M— — Ir sin2p+ I An Den gl +Dsin2ng+n-sin2(n+1) gl} (36). 


i 

Spannungsenergie: Der Spannungsansatz für o und das Querbiegemoment M gleichen 
den entsprechenden Funktionen für Biegung in der Krümmerebene nach (11) und (12), nur 
ist cos durch sin ersetzt. Dies bedeutet, daß die Spannungsenergie für die Normalspannungen 
den gleichen Aufbau wie (13) besitzt. Wenn wir hierzu noch die innere Arbeit der Gruppe II, 
die keine Freiwerte enthält, hinzuzählen, so erhalten wir die gesamte Spannungsenergie. Es 
genügt also, das Minimalprinzip auf die Spannungsenergie der Lastgruppe I anzuwenden. Aus 
der bereits festgestellten Übereinstimmung des Aufbaues mit der Spannungsenergie nach 
Gl. (13) erhalten wir für A, dieselben x-Werte, wie sie in Bild 6 dargestellt sind. Damit 
gilt für das Ringelement: 

M’cos’a-R-da M’sin’a:- R-da 

Erin 3 Am in 
M - cos a= Biegemoment quer zur Krümınerebene des Rohrquerschnitts, 
M - sin a = Torsionsmoment des Rohrquerschnitts. 
M’Rda, 
2Er’tn‘ 
Die Spannungsenergie eines Rohrbogens ergibt sich durch Integration über a. Die gegen- 
seitige Verdrehung der Endquerschnitte läßt sich nunmehr für jeden Belastungsfall aus der 
Spannungsenergie ermitteln. » wird in der gleichen Weise berücksichtigt wie unter B, II. 

Die Querverformung ergibt ein Oval, wie es bereits in Bild 2 dargestellt ist, jedoch 
sind die Symmetrieachsen um 45° gedreht. Infolge der Torsion tritt dann noch eine Ver- 
drillung hinzu, für die das Ringelement als gerader Zylinder von der Länge R-da aufgefaßt 
werden kann. 

Da die Spannungszustände der beiden untersuchten Fälle: Biegung in der Krümmer- 
ebene und quer hierzu, ähnlich sind, ist zu vermuten, daß auch hier der Einfluß von o=r/R 
auf die Verbiegung gering ist. Deshalb wollen wir auf eine genauere Untersuchung verzichten, 
da sie wesentlich umfangreicher ist, als die im Abschnitt B unter III durchgeführte. 507 


d A; = 


dA=dA,+dA,, = x».co®a+t(l+»)sina . . . .....+1(37). 
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Über die Analyse kurzer Kurvenzüge. 
Von Hans Blume in Hamburg. 


Es wird ein Verfahren zur Feststellung der Voraussetzungen dafür her- 
geleitet, daß eine kurze Registrierkurve durch eine Summe exponentiell 
abklingender Wellenzüge konstanter Frequenz dargestellt werden kann. 
Ein Weg zur Berechnung der dabei auftretenden Konstanten wird ge- 
geben. Da die Methode nicht auf die Darstellung durch solche Komponenten 
beschränkt ist, wird weiter die Frage erörtert, welche Arten von Wellen- 
zügen zweckmäßig einer Analyse zugrunde zu legen sind. 


Die Analyse von Registrierkurven nach persistenten und quasipersistenten Bestand- 
teilen ist mit Hilfe der Periodogrammrechnung möglich (Stumpff, Grundlagen und Methoden 
der Periodenforschung, Springer, Berlin 1937). Sie lehrt, eine reine Sinuswelle innerhalb 
einer Registrierkurve mit ihrer wahren Wellenlänge, Amplitude und Phase zu erkennen. 
Ihre beiden praktisch wichtigsten Hilfsmittel, das Amplituden- und Phasendiagramm, liefern 
aber nur dann richtige Ergebnisse, wenn der Geltungsbereich der zu analysierenden Kurve 
groß gegenüber der Länge des Analysenintervalls ist, und versagen in den Fällen, wo der 
Wellenzug „kurz“, d.h. kleiner als das Analysenintervall ist (Stumpff, l. c., S. 127). 

In folgendem wird ein Verfahren mit Hilfe der Periodogrammrechnung abgeleitet, 
durch das angegeben wird, unter welchen Voraussetzungen eine kurze Registrierkurve am 
zweckmäßigsten durch eine Überlagerung exponentiell abklingender Wellenzüge konstanter 
Frequenzen dargestellt wird. Ferner wird gezeigt, wie sich die charakteristischen Größen 
dieser Wellenzüge berechnen lassen. Die Brauchbarkeit der Methode wird an praktischen 
Beispielen demonstriert. Ihre Anwendbarkeit ist nicht auf die Heranziehung exponentiell 
abklingender Wellenzüge konstanter Frequenzen zur Darstellung komplizierterer kurzer 
Kurvenzüge beschränkt. Daher wird weiterhin gezeigt, unter welchen mathematischen Vor- 
aussetzungen durch das Analysierverfahren erkannt werden kann, welche Art Wellenzüge 
zweckmäßig zum Aufbau eines komplizierteren Kurvenzuges heranzuziehen ist. Das Ver- 
fahren ist auf die Analyse diskreter Beobachtungsreihen sinngemäß anwendbar. 


81. Die er eines kurzen Kurvenzuges in Kosinusform. 


Zunächst seien einige Definitionen gegeben, die im 








gs 2 IN /\ . folgenden benutzt werden: Es sei (Bild 1) f(f) der ana- 
a \/ V Y A Iytische Ausdruck für einen he Kurvenzug. Er 
oem ta» soll für einen positiven Bereich der Länge I definiert sein 
Bild ı. (0<st,stst„+l) und „kurz“ heißen, wenn / kleiner ist 


als das Analysenintervall p, mit dem die Analyse vorge- 
nommen wird. Ferner soll immer {„,+1i<p sein, d.h. das Analysenintervall p soll den Be- 
reich ! vollständig einschließen für jeden Anfangspunkt q des Analysenintervalls: 0<g<t.,. 
In Verallgemeinerung des Integrals für die Fourierkoeffizienten werden folgende Integrale 
für jeden Parameter x des Bereiches -—o<x< + gebildet: 


ta+l 
a(,=-, \re-Weost-dedt.....2.2.2..2.. de) 
ta 
ta+l 
&0=, \ BB Ba ii re rer. AR 
tu 
ta +1 
a +Hibd=, \ RN. ums au ı SEE | ° 
ba 
Setzt man für x die diskrete Zahlenfolge x = x,, - nm=1,2,3...) ein, so erhält man die ge- 


wöhnlichen Fourierkoeffizienten der Funktion f(t) für das Analysenintervall der Länge ». 
Unter dem Amplitudendiagramm des kurzen Kurvenzuges f(t) verstelit man den Ausdruck: 





KEERRBF FERN : . - ». :» 2: 2: 2:0. 


Hier denkt man sich A (x, g) für einen beliebigen festen Wert <=, als Funktion von g oder 
für einen festen Wert von q als Funktion von x dargestellt. Aus dem (Juotienten: 
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b(x,q) 
aa, srid a 3 


k(x, = 


bildet man das Phasendiagramm, indem man die zu y(x,g) gehörige Bogenlänge für ein 
festes <—=x, als Funktion von g darstellt. (Bei allen folgenden Untersuchungen ist der 
harmonische Analysator von Coradi benutzt. Da er den Cosinuswert mit umgekehrten 
Vorzeichen liefert, werden in den Diagrammen dieser Arbeit die theoretisch errechneten 
Neigungen von Geraden mit entgegengesetzten Gradienten wiedergegeben.) 


Nun sei: f(l)=ccos(at— ß). Nach Einsatz in (2) ergibt die Rechnung: 


sin (= — a) 1/2 


eillix—e)l2—xgq +@ta+ 2] 
(x — a) 1/2 


c 
a(z,g)+ib(xz,g)= ry I 
(2a). 
se sin (x + a) 1/2 


p (x + a) 1/2 


Nach Stumpff (l. e., S. 103) ist der zweite Summand als Störvektor vernachlässigbar klein, 


eilx+e)l2—xy—aty 3] 


2 
“N e . - . . 
wenn 123 z.; 50 daß unter dieser Voraussetzung das Amplituden- und Phasendiagramm allein 


vom ersten Summanden als Hauptvektor geliefert wird. Da in seinem Argument die Größe g allein 
mit dem Faktor — x auftritt, kommt also im Phasendiagramm die zugrunde liegende Frequenz a 
nicht zur Anzeige, wie es bekanntlich für !>p der Fall ist. Dies erreicht man jedoch auch 
für l< p, wenn man bei der Analyse folgendermaßen vorgeht: Es sei ,=p=1 und die 
Analyse werde für eine diskrete Wertefolge q=0, q,,9,... mit konstantem Analysenintervall » 
durchgeführt. Man schiebt also schrittweise das Analysenintervall p in den Bereich des 
Kurvenzuges hinein und führt nach jedem Schritt die Analyse mit dem harmonischen Ana- 
lysator durch. Dabei hat man drei Fälle zu unterscheiden: 


1. Das Intervall p dringt in den Bereich ! ein. Dieser Eintrittsvorgang hält an, solange 
0<g<i ist. 
2. Das Analysenintervall umschließt den Bereich ! des Kurvenzuges ganz. Dies gilt für 
alle !<gq<p. . 
3. Das Intervall p tritt mit seinem Ende in den Bereich } ein. Dies gilt für den Werte- 
bereich p<qg<sp+Hl. 
Für die drei Fälle muß der Verlauf des Amplituden- und Phasendiagramms betrachtet werden. 
Von ihnen ist der zweite Fall soeben behandelt worden. 
Für den Eintrittsvorgang ergibt sich in komplexer Schreibweise: 
p+4 p+qa 
\ eilix-a)t+pa-ya+3ldt + I \ eilix+e)t—pa—yax- Id 
e% 


J 


» D 


c 
alz,+ib(z,)= = 
c sin(@—a)g/2 ‘ 
in a vie! [-(x +e1gl2+px+?] 
p1 (2 — a) g]2 R a 

„2 mw +regi 

TI R+a)gl2 
Der zweite Summand kann als Störvektor vernachlässigt werden, wenn q so groß wird, daß 
es mindestens drei Wellenlängen der zugrunde liegenden Kosinuskurve umfaßt (Stumpff, 
l. ec. S. 103). Da der Eintrittsvorgang mit q=1 abschließt, soll vorausgesetzt werden, daß 





eil-ıx—-e)g 2+px-— 3] 


9 
1>3 —, so daß gegen Ende des Eintrittsvorganges der Einfluß des Störvektors vernach- 


lässigt werden kann. Das Phasendiagramm stellt im wesentlichen das Argument des Haupt- 
vektors von (5) in Abhängigkeit von q graphisch dar. In ihm kommt 4 nur linear mit dem Faktor 


1 j \ 27 
—(2,„+ a) vor: Für jede Harmonische x = = m —=1,2..) ergibt also das Phasendia- 


gramm eine Gerade mit der Neigung — (£„+a)>, so daß aus ihr und der bekannten Har- 


monischen x, die Frequenz a bestimmt werden kann. Beim Übergang in den Vorgang 2, in 
dem das Analysenintervall den gesamten Kosinuskurvenzug umschließt, muß an den Stellen 
q=1l ein Knick im Phasendiagramm auftreten, da für I<q=<p die Neigung bei jeder Har- 
monischen x,„ gleich ist —x,„. Dieser Knick ist um so stärker, je mehr die Versuchs- 
frequenz x,„ von der zugrunde liegenden Frequenz a abweicht, und er verschwindet, wenn 
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“„=a ist. Beim Übergang in den dritten Teil der Analyse erfolgt ein zweiter Knick im 


Plhasendiagramm, da beim Austrittsvorgang wieder die Neigung — (x„+ a) > auftritt. 


Aus denselben Gründen muß nach Gl. (5) für <=a das Amplitudendiagramm im Be- 


. ö i RT ö u 
reich O<q<I! von Null an bis zum Wert — ! steigen. Im Bereich I<q=p ist es konstant 
und fällt im Bereich p<q<p-+l in gleichem Maße ab, wie es im ersten Teil angestiegen 
ist. Aus den Knickstellen im Amplitudendiagramm ist ! bestimmbar und damit aus der kon- 


stanten Amplitudenhöhe ra die Amplitude c, da a aus dem Phasendiagraınm bekannt ist. 
9 


i R = ’ ü 2 ER a 5 Zn 
Die Phase $ ergibt sich, indem man die zu a gehörige Wellenlänge errechnet und nun 
S m ” a 


2n , En 
mit irgendeinem Intervall p = „nn= 1,2..) mit dem Anfangspunkt qg=t„=p die Kurve 


b (a, g) 


ergibt die gesuchte Phase nach Gl. (2a). 
a(a,dg) 


nochmals analysiert. Der Quotient k (a, g) = 


Dies Verfahren werde am Beispiel des kurzen Kosinuskurvenzuges der Länge != 240 mm 


> 
und der Form y=f(t) =50 sin (® ) mit einem Analysenintervall p =340 mm geprüft. Das 


Amplitudendiagramm ist in Bild 2 wiedergegeben. Man sieht, daß der Amplitudenverlauf 








Bild 3 (rechts). 








Bild 2 (unten). 
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der sechsten Harmonischen im konstanten Teil ein Maximum gegenüber dem der übrigen 
Harmonischen einnimmt, woraus nach Gl. (5) folgt, daß diese Versuchsfrequenz x,=6 der 
gesuchten Frequenz am nächsten kommt. Der Anstieg und Abfall ist leicht wellig infolge 
der Einflüsse des Störvektors. Die Länge des Anstiegs bzw. Abfalls beträgt 230 mm. Dies 
wäre also die aus dem Amplitudendiagramm abgelesene Länge des Kurvenzuges, welche in 
Wirklichkeit 240 mm beträgt. Das Maximum der Amplitudenhöhe ist 31,5 mm. Der Am- 
sin (x — a) q|2 
(2 — a) g|2 
beeinflußt, je weiter nach Gl. (5) die Harmonischen x, von der Frequenz a entfernt liegen. 
In den Amplitudendiagrammen der Harmonischen x,=5, und x,=4 kommt das deutlich zum 
Ausdruck. Im Phasendiagramm (Bild 3) ergibt sich der Wert 5,6, während die wahre 


plitudenverlauf der übrigen Harmonischen wird um so mehr durch den Faktor 


Frequenz «a mit 7x = 5,67 angezeigt wird. Die Berechnung der Amplitude ce geschieht nach 


60 
Bestimmung von ! und a aus der Gl. (5): c=51,5 mm. Auch hier liegt bis 3°/, Überein- 
stimmung vor. (Auf die Gründe für die Knickstellen im Phasendiagramm der dritten und 
vierten Harmonischen an den Stellen für q, an denen im zugehörigen Amplitudendiagramm 
ein Minimum liegt, soll nicht weiter eingegangen werden.) 
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$ 2. Die Analyse eines kurzen Kurvenzuges mit konstanter Frequenz 

und exponentiell abklingender Amplitude. 

Der Kurvenzug habe die analytische Form: y=f(t)=c e-?!cos(at— Pf). Die Analyse 
soll in derselben Weise vor sich gehen, wie im vorigen Paragraphen, d.h. das Analysen- 
intervall wird schrittweise durch den Geltungsbereich ! des Kurvenzuges hindurchgeschoben. 
Dabei enstehen wieder die drei Stadien der Analyse, nämlich das Eintritts-, Mittel- und Aus- 
trittsstadium mit denselben Geltungsbereichen für q. Für den Eintrittsvorgang lauten die 
verallgemeinerten Fourier koeffizienten folgendermaßen: 

Daraus ergibt sich in komplexer Schreibweise: 


af.r. q) +ib(e, 4) 


— * gq+ila(p—a)—eg+?] e'’[x(»-a 
"rs —o 


1 f 


, oeq+ilxip- D+agy-?] il (vd )— FIX 

- e-®rd l I I etii / / | 

5 x D 
o—+t(e ta) 





Den ersten Summanden nennt man wieder den Hauptvektor, den zweiten den Störvektor. 
Das Quadrat ihrer Absolutbeträge ist: 
e* 


va +1 —2e-27cosq(-a) . . .. . (6a), 


ey a fe 
pP’ [o + — a)?‘ 
hi= E (e-?va +1 —2e-tacosg(x + a)) (6b) 
’ p[o?+le + a)®] ‘ . : A a 7 x 
Zunächst wird festgestellt, unter welchen Bedingungen der Störvektor gegenüber dem Haupt- 
vektor vernachlässigt werden kann. Zu dem Zwecke wird in den Gl. (6a) und (6b) <=, 
> 
“N Oo . . .- 
1=Tr und e = . gesetzt und der Quotient gebildet. Man erhält: 


e* 1-+e-t7re -2e-?rrecoos(4 rr) 


ep — (be). 
(1 u e®—+4 j 


= 
Bei sehr kleinem Wert von e weicht der exponentiell abklingende Kurvenzug nur wenig von 
einer reinen Kosinuswelle ab, so daß dann praktisch der im vorigen Paragraphen unter- 
suchte Fall vorliegt. Für beliebige Werte von e==0 verschwindet d nicht mehr mit unend- 
lich werdenden Werten von r. Man muß dann für die Praxis der Analyse zwei Fälle 
unterscheiden: 


l. Der Wert e= - ist so klein, daß die Einwirkung des Störvektors auf den Hauptvektor 
( 


vernachlässigt werden kann. Für die Praxis dürfte es noch tragbar sein, wenn der 
(6) 

Absolutbetrag des ersteren höhstens 10°/, vom letzteren ausmacht. Dann muß e= -& 0,2 
( 

sein. Dieser Fall wird in diesem Paragraphen untersucht. 

Der Wert e>0,2. Dieser Fall wird im nächsten Paragraphen behandelt. 


Für den Geltungsbereich der zu untersuchenden Kurve würde es genügen, wie eine 
Auswertung von Gl. 6e lehrt, daß er mindestens drei Wellenlängen umfaßt, damit bei der 
Analyse am Ende des Eintrittsvorganges und während des Mittelstadiums, wo das Analysen- 
intervall den Geltungsbereich ! umschließt, der Störvektor praktisch vernachlässigt werden kann. 

Es bleibt noch zu prüfen, welche Einschränkungen für den zu untersuchenden, exponen- 


z R 2 0 ® 3 er i 
tiell abnehmenden Kurvenzug daraus folgen, daß e = > den Wert 0,2 nicht überschreiten 
C 


97: 
soll. Nun bedeutet e « den Faktor, um den die Amplitude verkleinert wird, wenn ? um 
den Betrag der Wellenlänge des kurzen Kurvenzuges wächst. Dieser Faktor ist für den 
angegebenen Wert « ungefähr 0,28. Zusammenfassend läßt sich also feststellen: Der Stör- 
vektor wird im Hiöchstfalle 10°, des Hauptvektors ausmachen, wenn die Amplitude des 
Kurvenzuges beim Fortschreiten um eine Wellenlänge nieht unter 28°/, ihres Betrages fällt 
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und der Kurvenzug mindestens drei Wellenlängen umfaßt. Diese Verhältnisse sind derart, 
daß ein erheblicher Bereich der Anwendungen von ihnen erfaßt wird. 

Unter diesen Voraussetzungen ergibt sich für das Amplitudendiagramm nach Gl. (6a) 
in der Nachbarschaft von qg<I!, wenn «=a ist, die Formel: 





[0 


c . 
ren aa Are a EEE EEE REFGEERE © 


Hieraus folgt, daß das Amplitudendiagramm für solche Harmonischen z=x, 
> 
= r n(n=1,2...), die in genügender Nähe von a liegen, sich in einem exponentiellen An- 


. . r c Fr n1 ® 4 .. 
stieg dem asymptotischen Wert De annähert, wobei zu beachten ist, daß diese Annäherung 


für q=1 praktisch erreicht ist. Da dieser Asymptotenwert in dem Amplitudendiagramm 
durch die Analyse erhalten wird, kann man in der Gl. (7) ihn einsetzen, wodurch man eine 
Gleichung in o allein erhält. Aus ihr ist o bestimmbar, indem man in dem Amplituden- 
diagramm zu bestimmten Werten von g die zugehörigen Werte abliest. Kennt man also die 
Frequenz a, so sind o und ce nach Gl. (7) berechenbar. 

Zur Bestimmung von a dient das Phasendiagramm. Allerdings kann dieses in dem 
Stadium des Eintrittsvorganges der Analyse nicht dazu herangezogen werden. Denn nach 
Gl. (6) hat hier der Hauptvektor unter den gemachten Voraussetzungen in der Nähe von 
q=1 die Form: 

. 1 eilp—)x+?3] 

—o+(x—a)i j 
In diesem Vektor kommt 4 linear mit dem Koeffizienten — x vor, d.h. das Phasendiagramm, 
welches im wesentlichen durch diesen Vektor wiedergegeben wird, ist für jede Harmonische «,, 
gegen Ende des Eintrittsvorgangs eine Gerade mit der Neigung — x„(n=1,2..). Die ge- 
suchte Frequenz « kommt also nicht zur Anzeige. 

Man kann beweisen, daß dies im Phasendiagramm des Austrittsvorgangs unter den 
schon gemachten Voraussetzungen für den kurzen Kurvenzug geschieht. Denn hier ist: 


i —# 
a(z,)+ibix,g)» > 





p+I 
u{ 
(a, +, | e-et-m+il«-at-zaten+2ldt 
a 
p+I 
c . 
+ — \ e-rtt— p+ila+a)t—xqa—ap—3]ldt ! 
a (8a). 
q 
c 1 . _ 
FE A vw i 7) ]— e-e(a-p)+il-eg+ter+ J} 
»?—o —a 
c 1 ; j ; 
A SEE 7, ] — e-.ı—p)+ilegq—ep+ } 
p-o+tlc+ta 


Der Austrittsvorgang beginnt mit g=p und umfaßt den Bereich p<q<p-+l. Unter den 
gemachten Voraussetzungen können der Störvektor und der erste Summand in der Klammer 
des Hauptvektors vernachlässigt werden, so daß für das Phasendiagramm in der Umgebung 
vom 9 >p im wesentlichen folgender Vektor bestimmend ist: 

—ce-1!-») 


pl-o+(x —a)i] 


Ein vollständiger Beweis hierfür wird in folgendem Paragraphen erbracht. In Gl. (8b) er- 


alz, +ib(a,g= erereHPfl ,. . . .. ., 


6) 
scheint für jede Harmonische «= x, = zn im Argument die Größe y linear mit dem Ko- 


effizienten —a, d.h. in der Umgebung von q>p entsteht für jede Harmonische eine Gerade 
mit der Neigung — a, und damit ist a bestimmbar geworden. Die Anzeige von « tritt noch dadurch 
besonders auffällig in Erscheinung, daß das Phasendiagramm im allgemeinen an der Stelle y=p 
einen Knick besitzt, der an der Stelle 4=1 nicht auftritt. Denn wie beim Phasendiagramm 
im Bereich des Eintrittsvorganges bewiesen wurde, ergibt sich gegen Ende des Bereiches 0O<q<I 
für jede Harmonische ©, ein Gradenstück mit der Neigung —x,„. Sobald der Eintrittsvor- 
gang beendet ist, geht die Analyse in das Stadium über, in dem der Geltungsbereich ! des 
kurzen Kurvenzuges ganz umschlossen ist. Dann gilt für den verallgemeinerten Fourier- 
koeffizienten in komplexer Form: 
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p+I 


Sr c ; a i + 
az ril,d=, \ a !—-p)+ifar— e)t—xga+ep dt 
) 


» 


p-+i 
c —o(t—p)+ilx +eit— xa—ep—?}] 
+— | emet-ı Zara. dt 
Pr”, (9). 
v 
o 1 l 1 l 3 ap +ß 
u EN : fe- o+ijix-a)\(p+D—xqa+ap+Pß] _ eil(ix—a)p—xqtap+ In 
» —-grria—y 
= 1 





Pat nen (et +ilea+ e\(p+ND—xq-—ep—?| __ eilix +a)p— xq—ap— 31 
p-o+tilc+ta) , 
Unter den gemachten Voraussetzungen ist hier ausschlaggebend allein der Vektor: 

1 


-— N. us sack VRREEE n 
-P +i(x — a) 


—c 
alz,)+tib(x,g) » > 
Seine Amplitude ist von q unabhängig und infolgedessen das Amplitudendiagramm jeder 
Harmonischen x„ konstant. Das Argument ist linear in g mit dem Keoffizienten — x, für 
jede Harmonische x&„. Daher kommt im Phasendiagramm die Frequenz «a nicht zur Anzeige 
und an der Stelle q=] tritt kein Knick auf. Da aber von g=p an die’ zugrunde liegende 
Frequenz a angezeigt wird, muß an dieser Stelle im Phasendiagramm ein Knick auftreten, 
Er wird um so stärker sein, je weiter die Versuchsfrequenz x„ von a abweicht, und ver- 
schwindet für z«„=a. Diese Verhältnisse lassen eine genaue Bestimmung von a aus. dem 
Phasendiagramm im Stadium des Austrittvorganges der Analyse erwarten. 
Da mit der Kenntnis von a die Berechnung von o und c möglich ist, wie schon gezeigt 
wurde, so sind die charakteristischen Größen a,o,c, ö des kurzen Kurvenzuges bestimmbar. 
Das Verfahren wird an einem Beispiel erläutert und geprüft. Es sei (Bild 4): 
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Hier ist e&= —-= 0,22. Da diese Kurve mit p = 168mm analysiert wird, so muß die Frequenz a 


a 
mit der Größe 168:34=4,94 im Phasendiagramm zur Anzeige kommen. Bild 5 zeigt das 
Phasendiagramm. Hier kommt der theoretisch abgeleitete Verlauf, insbesondere die Knick- 
stelle und die Parallelität, klar zum Ausdruck. Für die Frequenz ergibt sich aus der Neigung 
der parallelen Stücke der Wert 4,8. Im Amplitudendiagramm (Bild 6) ist der besseren Über- 
sicht wegen nur diejenige Harmonische x,=5 gezeichnet, welche der wahren Frequenz 
am nächsten liegt und daher den höchsten Asymptotenwert von 7,25 mm erreicht, d. h. 
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Nach Gl. (7) erhält man für den Amplitudenverlauf im Stadium des Eintritts- 


pP oO 
vorganges für o die Gleichung: 


1,25. 
a EB ed 2} GE 4? 


Die bei der Analyse auftretenden Schwankungen im Amplitudenverlauf sind auf die Ein- 
wirkung des Störvektors zurückzuführen. 
Die Bestimmung des Exponentialkoeffizienten mit Hilfe der Gl. (7) aus dem Amplituden- 
verlauf ergibt folgende Wertetabelle: 
q 23132161904 83256106416 
0-10”? 431 38137/14144 | 44 | 40 | 42 |42 | 38 | 39 





Als Mittelwert ergibt sich o=0,042. Man sieht, daß man einen ziemlich genauen Wert für o 
erhält, ohne daß man sich auf einen zu engen Bereich für 4 beschränken braucht. Durch 
N. » c - = fe} a - y- - 
Einsetzen von o in — =7,25 errechnet man: ce = 168 - 0,042: 7,25 = 51 mm. 
po 
Zusammenfassend läßt sich feststellen, daß durch das abgeleitete Analysierverfahren 
die charakteristischen Größen der zugrunde liegenden Kurve mit einer Genauigkeit von 2°, 


bestimmt werden können. 


$ 3. Die Trennung zweier verschiedener kurzen Kurvenzüge 
mit exponentiell abklingender Amplitude. 
Die analytische Form der Kurvenzuges sei: 
y=ftll) =c, e-tıteos(a,t— B)+w,e-ttcos,t—P.) - » -» :.. MM. 


Der verallgemeinerte Fourierkoeffizient dieser Summe ist gleich der Summe der verallge- 
meinerten Fourierkoeffizienten der beiden Summanden. Setzt man zunächst voraus, daß 
für jeden Summanden die [im $ 2] gemachten Voraussetzungen gelten, so wird für den Aus- 
trittsvorgang in der Analyse der verallgemeinerte Fourierkoeffizient im Bereich p<q<p+HI 

folgenden Wert besitzen: 
i —c, e=Pılqa PP j u 2 
al, « l ib x,g)» —; eil-ıa+ up+Pı] 

ET 0, til —a,) 
—02(9 —») 

€, n eil @4 +02» + 32] 

p 0, +?(x Q,) 

Zur Abkürzung sei gesetzt: 

1 1 7 (x — a,) 
;‚e®, wo tgy= : 


- =o- ; j ist, 
"H»trlz—a,) yo,rl&—a,) 0» 


(12a). 
6, e-%(q—p») 


> W=E-,4t,Pp+h tm M=WwtRT, (r 


yvo,+(x — a,)' 


| 


Dann erhält man: 
az, ya +iba,g-h,kef'rn.e". . . 2. 2 2 222. (12b). 
Ohne Einschränkung der Allgemeinheit wird angenommen,’ die Ver- 
teilung der Indizes sei so gewällt, daß h,>h,. Die Vektor- 
summe (12b) wird in einer komplexen Ebene aufgetragen (Bild 7). 
Durch Anwendung des Sinussatzes erhält man nach Stumpff, 1. e., 
S. 122: 
h,sin(p,— 9,) Zr 
NE TE I 
h,+h,cos(y, —9,) 


tgö= 


Es sei zunächst folgendes angenommen: Wenn 4 von p bis p-+I 
wächst, so bleibe in diesem Bereich des Austrittsstadiums der 
Analyse h, >h, für dieselbe Frequenz. In der Vektorebene dreht 
sich dann der Vektor h, ei: mit der Winkelgeschwindigkeit — «a, 
um den Nullpunkt, während der Vektor h,eif: mit der Wirkel- 
geschwindigkeit — a, um den Endpunkt des ersteren rotiert. Der resultierende Vektor bildet 
mit dem ersteren den Winkel ö. Dieser Winkel besitzt nach Gl. (13a) ein Maximum für 
den Fall: 





Bild 7. 
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Der Wert dieses Maximums ist nach Gleichung (13a): 
h, 


eo = = 
Sl Sud Tea 1 


Dieser Ausdruck ist reell, da h, >h,. Der Winkel ö,„ ist immer kleiner als z/2 und hat im 
Grenzfall k,=h, den Grenzwert z/2 für 9,=g9,+7. Daraus folgt: Der Winkel des resul- 
tierenden Vektors mit der waagerechten Achse in der Vektorebene ist y()=y,(Q)+ö(g) 
im Bereich p<gqg<p-+l. Dieser Winkel besteht also aus einer Überlagerung des Winkels 
ö(g) über den Winkel @,, der sich seinerseits mit der Winkelgeschindigkeit — a, gleichförmig 
um den Nullpunkt dreht. Trägt man daher die Bogenlänge von y(g) als Funktion von q auf, 
wie dies im Phasendiagramm geschieht, so erhält man eine Wellenlinie, deren Mittellinie mit 
der g-Achse den konstanten Winkel —a, bildet und deren Periodenlänge @ beträgt: 
an _ 2 
d, —u,  2n 
_ (n,—n;) 
PP 

Bild 8 zeigt schematisch den beschriebenen Verlauf des Phasen- 
diagramms an. Ist die Periodenlänge @ mindestens einmal ganz 
im Bereich p<g<p-+l enthalten, so ist also a, aus dem 
Phasendiagramm grundsätzlich bestimmbar. Dies wird praktisch 
um so besser möglich sein, je kleiner h, gegenüber h, ist. 

Hieraus folgert für die Analyse das wichtige Ergebnis: 

Liegt der Analyse ein kurzer Kurvenzug mit dem analyti- 
schen Ausdruck der Gl. (11) zugrunde, so wird im Austritts- 

>} 
h A 2n a 
stadium des Phasendiagramms jede Harmonische «= x, = r n(n=1,2,...) diejenige von 
) 


den beiden zugrunde liegenden Frequenzen a,, a, anzeigen, für welche die Amplitude A, in 
der Vektorebene die größere ist: 


C, e-%,(4- pP 


PV» tler —a,) 


Dies besagt zunächst einmal, daß beim Vorliegen eines einzigen Kurvenzuges nach $ 2 
der Störvektor die Frequenzbestimmung aus dem Phasendiagramm des Hauptvektors grund- 
sätzlich nicht beeinträchtigen kann, wie groß auch e=o:a sein mag. Denn in diesem Falle 
ist = — a, &,=C,, 0%,=o, und infolgedessen der Wert h, des Hauptvektors immer größer 
als derjenige des Störvektors, der nur positive Werte.besitzt beim harmonischen Analysator. 

) 
F ‚ F 7 Be . " j i 
Außerdem ist der Austrittsvorgang von der Länge 1>3—=6—-, während die Periode der 
a, a, 
im Phasendiagramm auf der Harmonischen x„=a, erscheinenden Störwelle die Länge 

9 

=, 7T r r . . . 
0=- we besitzt. @ ist also achtmal im Bereich des Austrittsvorganges enthalten. 

u a, 


' 


v=1,2). 


i x 1 
Sind darüber hinaus noch die Voraussetzungen von $ 2 über & erfüllt, so ist h,< 10 h, und 


infolgedessen die Anzeige der Neigung a, im Phasendiagramm praktisch ungestört. Von dieser 
Tatsache wurde im Beispiel am Schluß des letzten Paragraphen Gebrauch gemacht. 


Ist für eine Harmonische die Voraussetzung h,<h, nicht erfüllt, sondern h,=h,, so 
kommt im Phasendiagramm dieser Harmonischen keine der briden Frequenzen zum Vorschein. 
Sie trennt die übrigen Harmonischen x, in zwei Klassen, von denen die Angehörigen der 
einen die Frequenz a, und die der andern a, anzeigen, so daß dennoch grundsätzlich beide 
Frequenzen bestimmt werden können, wenn genügend Harmonische vorliegen. Dies wird ins- 
besondere dadurch immer erreicht werden können, daß man das Analysenintervall genügend 

9 


groß wählt, damit die Werte vieler Harmonischen «= r,= > n(n=1,2...) zwischen den 
5 2) 7 
n, und .=7 n, liegen. Ist z. B. c,=ct,, 09, =0,=0, so wird: 


Frequenzen a, = 


h,_ı/(po®+4n’(n—n,)” 
h, Y (po®+4n?(n—n,)? nu Fuath de 


" u n,+n, A i 4 ’ 
Für alle Harmonischen zu (n >-5 ") wird die Frequenz a, bzw. n, und für die Harmo- 


24 
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n, + 


\ n i : ' ' 
nischen zu [m < HZ) die Frequenz a, bzw. n, angezeigt. Durch die Harmonische 


En (n= ud wird weder a, noch a, erkennbar sein. 

Zusammenfassend kann man feststellen, daß es in der Praxis im allgemeinen möglich 
sein wird, die beiden zugrunde liegenden Frequenzen ziemlich genau zu bestimmen, wenn man 
einen Analysator mit hoher meßtechnischer Präzision verwendet, wie dies z. B. bei dem 
Analysator von Coradi der Fall ist. 

Die Ableitungen über die prinzipielle Trennbarkeit zweier Frequenzen werden ver- 
wickelter, wenn im Bereich des Austrittsvorganges die Größenverhältnisse h, und %#, für eine 
bestimmte Harmonische x, nicht konstant sind, sondern wechseln. Dies zwingt in der Praxis 
der Periodogrammanalyse zur Vorsicht, beeinträchtigt sie jedoch, wie die Erfahrung zeigt, in 
keiner Weise entscheidend. 

Hat man die Frequenzen a, und a, auf die angegebene Weise bestimmt, so kann man 
durch passende Wahl des Analysenintervalls grundsätzlich erreichen, daß jede der beiden 
Frequenzen mit einer Harmonischen zusammenfällt. Liegen diese Frequenzen genügend weit 
auseinander, so bestimmt man aus dem Amplitudenverlauf dieser Harmonischen die Werte 
o, und c„ (r=1,2) in der Weise, wie es im zweiten Paragraphen dargestellt wurde, ohne 
auf die gegenseitige Beeinflussung der beiden Wellen Rücksicht zu nehmen. Sind beide 
Wellen benachbart, so berechnet ınan zunächst die Größe o, wie im zweiten Paragraphen 
angegeben ist, da in diesem Falle meistens, z. B. bei Stoßvorgängen, zu erwarten sein wird, 
daß diese beiden Exponenten im wesentlichen übereinstimmen. Diesen gemeinsamen Wert 
legt man als Exponenten für beide Frequenzen der weiteren Rechnung zugrunde. Bei der 
nun folgenden Bestimmung der beiden Amplituden c, hat man jedoch die gegenseitige Be- 
einflussung der Wellen mit den Frequenzen a, und a, zu berücksichtigen. Man geht dabei 
wieder von Gl. (12b) aus. Für den resultierenden Vektor ergibt sich nach Bild 7 durch An- 
wendung des Kosinussatzes: 

”(z,)=h(&,g)+h;(z,g)+2h, (x, g)h,(x, q) cos p' V=9,—Pp, (14). 
Da a,, a,, o, und o, bekannt sind, ist @’ für jeden Punkt des Phasendiagramms berechenbar. 
Besonders einfach und daher für die Praxis wichtig ist es, Stellen g im Phasendiagramm 
zu finden, an denen ’ gleich einem ganzzahligen Vielfachen von x ist, weil dann die rechte 
Seite von Gl. (14) ein reines Quadrat wird und somit eine lineare Gleichung vorliegt. Dieser 
Fall wird an dem Beispiel der Analyse einer Stoßkurve im folgenden Paragraphen dargelegt. 
Grundsätzlich erreicht man aus Gl]. (14) zwei Gleichungen mit den beiden Unbekannten c, 
und c,, indem man der Reihe nach x,=a, bzw. a, setzt. Dabei ändert sich cos’ nicht. 
Allerdings ist vorausgesetzt, daß man die Phasen £, und £, der Originalkurven nach GI. (11) 
kennt. Dies ist bei Stoßkurven der Fall, da hier erfahrungsgemäß ß, = ß, = W°* ist. 

Bisher war in diesem Paragraphen vorausgesetzt, daß der Einfluß des Störvektors auf 
den Hauptvektor vernachlässigbar klein sei. Ist das nicht der Fall, so kann es sein, daß die 
beiden Frequenzen a, und a, so weit auseinanderliegen, daß sie sich nicht wesentlich stören. 
Dann können mit den aus den Phasendiagrammen bestimmten Frequenzen die Größen o, 
und o, berechnet werden, wie im zweiten Paragraphen angegeben ist. Bei der Berechnung 
von c, und c, berücksichtigt man dann nicht die gegenseitige Beeinflussung, sondern nur den 
des Störvektors unter Anwendung von Gl. (14). 

Verwickelter werden die Verhältnisse, wenn weder der Einfluß des Störvektors jeder 
Frequenz noch ihre gegenseitige Beeinflussung vernachlässigt werden können. Die Frequenz- 
bestimmung ist jedoch auch in diesem Falle genau durchführbar, wenn für beide Frequenzen 
nach Gl. (9a) der erste Summand im Haupt- und Störvektor vernachlässigt werden kann, 
d. h. wenn e-!? genügend klein ist. Dies dürfte für einen erheblichen Bereich der Praxis 
erfüllt sein. Die Bestimmung von o, c, und c, wird ungenauer, da ja dem reinen Amplituden- 
abfall nach dem Exponentialgesetz sich mehrere störende Periodizitäten überlagern. In der 
Praxis braucht man hierbei jedoch nicht so ängstlich zu sein, wie das auf den ersten Blick 
scheinen könnte, wie das Beispiel im nächsten Paragraphen zeigen wird. 


8 4. Die Analyse eines experimentell gemessenen Stoßvorgangs. 


Die Verwendbarkeit der dargelegten Anaiysiermethode soll an einer Registrierkurve, 
die durch die Messung eines kurzen Stoßvorganges mit einem Schwingwegmesser erhalten 
wurde, erprobt werden. Die Kurve') stellt die vertikale Schwingung einer Betonstraße dar, 
welche dadurch entsteht, daß ein Lastwagen mit Vollgummireifen über einen fünf Zentimeter 
hohen Keil fährt. Die Originalkurve ist in Bild 9 wiedergegeben. 


1) Z. Instramentenkde. Bd. 62 (1942), S. 361. 
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14). Für die Analyse wurde sie durch optische Abbildung auf einen Bereich von 1 = 143 mm 
ar. vergrößert und mit einer Analysenintervallänge von p=177 mm analysiert. Bild 11 gibt das 
nm dabei erhaltene Phasendiagramm der vier Harmonischen #,=3, #,=4, 2,=5, x,=6 wieder. 
hte Auf der waagerechten Achse sind diejenigen Werte q verzeichnet, die bei der schrittweisen 
ser Analyse der Reihe nach als Anfangspunkte des Analysenintervalls p genommen wurden. Im 
gt. Bereich von O0 bis 190 mm hatten diese Punkte einen Abstand von 10 mm, von da bis zum 
c, Schluß der Analyse einen solchen von 5 mm. Als Ordinaten sind die Bogenlängen des Winkels 
| > 
- y(n, q) aufgetragen, wobei tgy(n,gq) = (m 9) für jede Harmonische <= x, = "nln=1,2.. .) 
) a(n,dg) p 
ist. Der Wert b(n,g) ist der am Analysator abgelesene Sinuswert für den Fall, wo das 
ruf Analysenintervall von p=177 mm im Punkt q beginnt, a(n,g) ist der zugehörige Kosinus- 
lie wert. Die Tatsache, daß im Phasendiagramm für alle Harmonischen gradlinige Fortführung 
en. über g=180mm bis 9=19% mm eintritt, ist darauf zurückzuführen, daß zu Beginn der 
0, Registrierkurve der Einschwingvorgang wirksam ist. Während jedoch von q9=19 an diese 
ng Gerade bei der vierten und sechsten Harmonischen verlassen wird, setzt sie sich bei der 
en dritten bis zum Wert q9=210 und bei der fünften bis g=220 fort. Bei der dritten Harmo- 
nischen ergibt eine rechnerische Neigungsermittelung zwischen 225 und 300 den Wert 2,97, 
er d. h. nach einer Phasenverschiebung zwischen 210 und 225 setzt sich die frühere Neigung 
1Z- fort. Der unregelmäßige Verlauf von =300 an bei allen Harmonischen beruht darauf, daß 
en die Meßgenauigkeit des Analysators unterschritten ist. Die dritte Harmonische besitzt das 
n, größte Maximum. Auf Grund der dargelegten Theorie wird angenommen, daß die durch sie 
iS angezeigte Frequenz n—3 im gesamten Kurvenzug wirklich vorkommt. Sie ist gleich 13 Hertz. 
n- Der Amplitudenverlauf Bild 10a und 1Ob der fünften Harmonischen im zweiten Stadium der 
er Analyse 160 <9< 1% zeigt das zweitgrößte Maximum, während im gleichen Bereich die da- 
°K zwischenliegende vierte Harmonische das kleinste Maximum besitzt. Da die Neigung der 
fünften Harmonischen sich bis 4=220 fortsetzt, wird angenommen, daß auch die ihr zu- 
gehörige Frequenz von 22,5 Hertz in der Kurve wirklich vorkommt. Im Bereich 220 < 9 < 300 
zeigt diese Harmonische eine andere Neigung an, deren rechnerisch ermittelter Wert 4,1 be- 
e, trägt. Im gleichen Bereich des Amplitudendiagramms erkennt man, daß die Amplitude der 
nn vierten Harmonischen die der fünften erheblich übersteigt. Daraus kann geschlossen werden, 
r, daß die Frequenz von 22,5 Hertz nur noch schwach im Auslauf der Kurve vorkommt und 
r ihre Anzeige von der stark auftretenden Frequenz 18,5 Hertz, die dem Neigungswert 4,1 ent- 





sprielit, unterdrückt wird, wie früher theoretisch begründet wurde. Das Phasendiagramm der 
vierten Harmonischen verläßt in krassem Sprung die Neigung n = #4 und zeigt die Tendenz, 
im Bereich 200 <9<220 sich der Neigung n —=5 anzupassen und am Ende dieses Bereiches 
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in eine Neigung n=3,8 (16,5 Hertz) überzugehen. Tatsächlich liegt ihre Amplitude in diesem Der 
Bereich tief unter allen anderen. Nach der Theorie muß man daher annehmen, daß die so ( 

Frequenz 16,5 Hertz vernachlässigbar klein im Kopf der Kurve vorkommt. Der kurze über- 

große Anstieg der vierten Harmonischen wird dabei unberücksichtigt gelassen. Demgegenüber 

kommt offenbar die Frequenz 16,5 Hertz im Auslauf der Kurve vor, wie schon besprochen 
wurde. Das Phasendiagramm der sechsten Harmonischen ist ab q= 19 unregelmäßig, ihre Zun 
Amplitude fällt schnell zu kleinsten Werten ab. Daher wird angenommen, daß keine wirk- tion 
lich vorkommende Frequenz von nennenswerter Amplitude zu ihr gehört. Zusammenfassend geg 
wird festgestellt, daß in dem hauptsächlich interessierenden Teil, dem Kopf der Kurve, die tigs 
beiden Frequenzen 13 und 22,5 Hertz bestimmend sind. Diese sollen daher allein berück- stüc 
sichtigt werden. ist 
Die nächste Aufgabe ist die Bestimmung des Exponentialabfalls für diese beiden Fre- man 
quenzen. Zu dem Zweck zieht man die Gl. (8) heran, die den Anstieg der Amplitude im Ein- Max 
trittsstadium der Analyse wiedergibt. Für die dritte Harmonische ist der asymptotisch er- sieh 
reichte Höchstwert h„;= 11,80, so daß zur Bestimmung von o, die Gleichung vorliegt: sch\ 
ad, )+ibßgl=118 (l—e-te). eine 
Ihre linke Seite kann aus dem Amplitudendiagramm zu jedem Wert q abgelesen werden. - ( 
Die Errechnung von o, ergibt folgende Wertetabelle: ge 
= » „ } ’ ni Ste 

q 20 | 30 | 40 | 50 | 60 | 70 80 1%  100| 110 120 130! 140 | 150 
0,10”? | 1,96 3,15 2,93 3,54 | 2,96 | 2,24 | 1,85 | 1,85 1,70 | 2,56 | 2,02 | 2,16 3,30) — 

Der Mittelwert ist 0, = 0,025. Auf dieselbe Weise verfährt man mit der fünften Harmonischen, 

} 6, ; 2 : . i 

wobei als Asymptotenwert 7 —10,1 angenommen wird. Dies ist der’ Wert, den die 

9; 

Amplitude im zweiten Stadium der Analyse besitzt. wo das Analysenintervall den gesamten 

Kurvenbereich umfaßt. Die im Anfangsstadium der Analyse erhaltenen Werte, welche größer 
als dieser Wert sind, scheiden bei der Bestimmung von o, aus. Dann erhält man folgende vorg 
Tabelle: kurı 
q | 20 | 30 | 40 | 50 | 60 | 70 | 80 | 90 | 100 ii 
0,:10-? 2,08 | 2,62 | 3,33 | 4,46 | 3,26 | 2,92 | 2,82 | 2,95 | 2,92 ' Ges: 
Fur Bye | die ] 
Aus ihr ergibt sich 0,=0,030. Der weiteren Rechnung wird für beide Frequenzen derselbe gese 


ar u IN - > leo 
Wert go =0,028 zugrunde gelegt. mat 


Zur Bestimmung der beiden Amplituden c, und ce, dient die Gl. (14), wobei es zwecks und 
Vereinfachung der Rechnung wichtig ist, eine Stelle zu finden, an welcher der Wert von 
cosp’'—=1 ist. Zunächst wird angenommen, daß die Phasenlage jeder der beiden abklingenden 
Wellen gleich 90 Grad ist. Dann ist: Her: 
=, -9M,=-@,-a)p- ty, —y))- ii 
Der Einsatz der bekannten Werte ergibt y,=0 und y,—=—-110°. Daraus folgt, daß nach Dee 
einer Phasenverschiebung — (a, —a,)(p—gQ)=110° der Winkel 9'=0 ist. Dies entspricht Ken 
einem Wert p—q von —27. Es wurde schon erwähnt, daß der Einschwingvorgang die ge- 
naue Größe des Wertes q, für den im Sinne der dargelegten Theorie der Austrittsvorgang 
beginnt, nicht bekannt ist. Aus dem Verlauf der Registrierkurve kann man entnehmen, d 
Fr 2 - in n \ er 
daß dies zwischen q gleich 180 und 185 der Fall ist. Diese Annahme wird durch andere u 


Überlegungen noch bestätigt werden. Hieraus folgt, daß der Wert q, für den die Phasen- Ba 
verschiebung 9’=0 beträgt, bei beiden Amplitudendiagrammen »=3 und n=5 in der Nähe 
von 49=210 liegen muß. Tatsächlich besitzen beide Diagramme hier einen Wendepunkt. 
Diese Tatsache kann als eine Bestätigung dafür angesehen werden, die gesuchte Stelle richtig 
gefunden zu haben, und zwar aus folgendem Grunde: Nach Gl. (14) besteht das Quadrat des 
Amplitudendiagramms aus der Überlagerung eines aperiodischen Teiles A? +h? und einem Bild 
Teil 2h, h,cos p', der mit der Periode von @’ exponentiell abklingt. In letzterem ist die er 
Differenz für zwei sehr benachbarte Werte von 4 an den Stellen seiner Maxima bzw. Minima am 
kleinsten. Diese liegen aber sehr nahe an den Stellen, an denen cos 9'= + 1 ist. Andererseits kann 
aus dieser Übereinstimmung, die sich aus der Rechnung unter Verwendung der ermittelten 
Werte für «,, a, und o ergeben hat, geschlossen werden, daß man diese Größen annähernd 
richtig bestimmt hat. Unter diesen Bedingungen, d. h. für py— q=--?27, erhält man aus 
Gl. (14) folgende beiden linearen Gleichungen zur Bestimmung von ce, und e,: 
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Der Einsatz der Werte für 9g=205 aus dem Amplitudendiagramm liefert e,=385; c,=35, 
so daß die Analyse folgendes Ergebnis hat: 


t) + 35 e-&08t sin 2n nt) ee an. 
177 


/ 


y=38,5e-"®t sin (2 a 


177 
Zum Vergleich mit der Originalkurve ist die Funk- ylımm) 
tion (15) im gleichen Maßstab in Bild 12 wieder- 3r 

gegeben. Vorzüglich ist die Wiedergabe des wich- #- 

tigsten Teiles der Kurve, nämlich des Anfangs- x» 
stückes bis zum ersten Nulldurchgang. Im übrigen 
ist der Charakter der Kurve wesentlich erfaßt, wie % 





man aus einem Vergleich der Lage gleichartiger v 
Maxima, Nullstellen und Minima beider Kurven er- ar 


sieht. Zu dem Zwecke muß man von dem Ein- Bild 12. 
schwingungsvorgang der Originalkurve absehen und 

einen andern Anfangspunkt bestimmen. Deswegen trägt man den Wert der ersten Nullstelle 
in Gl. (15) von der ersten Nullstelle der Originalkurve nach links ab. Man erhält so den 
Punkt O als Nullpunkt in der letzteren. Die Aufeinanderfolge der genannten charakteristischen 
Stellen in beiden Kurven wird durch folgende Tabelle dargestellt: 


Synthetische Kurve | O0 | 12 | 22 | 30 | 45 | 56 | 67 | 77 | 86 |105 
Originalkurve 0/10/2,30|43|57)67 | 77|8|10 





$ 5. Schlußfolgerungen und Erweiterungen. 


Aus den vorhergehenden Darlegungen ergibt sich für die Praxis der Analyse von Stoß- 
vorgängen, daß es im allgemeinen nicht gerechtfertigt ist, den ersten Hauptbogen der Registrier: 
kurve als sinusförmige Annäherung einer bestimmten Frequenz aufzufassen und diese dann 
als sogenannte „Stoßfrequenz* bzw. „Stoßamplitude* für die Umrechnung auf Schwingweg, 
Geschwindigkeit und Beschleunigung zu benutzen. Die bisher festgestellte Unmöglichkeit, 

> oO 
die Registrierkurven eines Stoßvorganges, der gleichzeitig mit einem Schwingweg-, Schwingungs- 
” . . 0 ° g . e . ° > 
geschwindigkeits- und Schwingungsbeschleunigungsmesser aufgenommen wurde, in eine mathe- 
matische Zuordnung bringen zu können ’?), dürfte mit auf diesen Umstand zurückzuführen sein 
und ist Gegenstand weiterer Untersuchungen. 


Das formale Prinzip der geschilderten Analysiermethode ist nicht auf die hier behandelte 
Heranziehung exponentiell abklingender Wellenzüge konstanter Frequenzen als Aufbauelemente 
für komplizierter abklingende Kurvenzüge beschränkt, wenn diese auch ihren praktisch 
wichtigsten Fall darstellt. Die Anwendung auf linear abfallende Wellenzüge konstanter 
Frequenz zum Aufbau komplizierterer Kurvenzüge verläuft analog. Das zugrunde liegende 
formale Prinzip läßt sich folgendermaßen ausdrücken: Es sei 

EEE net er A 
der Typ der Funktionen, die zum Aufbau eines komplizierteren, kurzen Kurvenzuges dienen 
sollen. Da diese Aufbaufunktion einfach sein soll, ist es hinreichend, von der reellen 
Funktion f(£) vorauszusetzen, daß sie für positive wachsende Werte von & eine stetig ab- 
nehmende Funktion sei. Ihr Integral wird dann bei Variabilität der einen Integrationsgrenze ' 
und Konstanthaltung der anderen eine Funktion folgender Form sein: 


y=F(£&,y .eivst®), 


Bildet man mit (16) in Analogie zu Gl. (8a) für das Austrittsstadium der Analyse den ver- 
allgemeinerten Fourierkoeffizienten, so erhält man: 


p+i 

1 ’ 
al, g)+ib(z, un \ ftt—p) dla -at-zarentPlgg 
i . (16a). 
p+l 


+, \ ft-peilata-xe-er-dlgı 
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Das erste Integral läßt sich auf folgende Form bringen: 
v 7 ı 


u; ri a 
N @ep+ß] \ ftp) e*- )edt | 
q (16b). 
1 it- +ap+ß ’ i —ea)(p+D+k ix —e)qg+k 
eh xq p+ Fl, x — a): eil‘“ )(p 4 I Fig —P,2 — a): eil‘“ 1 kN 
Ist Fl) genügend klein, so kann man für (16b) folgenden Ausdruck setzen: 


a). eil-«a tap+?+ kl 


Hier tritt im Argument die Variable q linear mit dem Faktor —«a auf, dies ist aber die 
Voraussetzung dafür, daß im Phasendiagramm für jede Versuchsfrequenz x die wahre 
Frequenz a als konstante Neigung zur Anzeige kommt. Eine hinreichende Bedingung dafür, 
daß F(l) genügend klein wird, liegt vor, wenn f(£) in Gl. (16) mit wachsenden positiven 
Werten von £, in genügend hoher Ordnung Null wird. Übrig bleibt dann noch eine Größen- 
abschätzung des Störvektors, d. h. des zweiten Integrals in Gl. (16a), das von derselben Form 
ist wie Gl. (16b), da lediglich die Größe a mit umgekehrten. Vorzeichen zu nehmen ist. 
Diese Überlegungen zeigen, daß die beschriebene Analysiermethode grundsätzlich er- 
weitert werden kann, und daß man aus dem Amplitudenverlauf der Harmonischen x =x,„=a 
grundsätzlich diejenige Funktion f (£), welche als Aufbauteil in der zu analysierenden Registrier- 
kurve steckt, unter verhältnismäßig allgemeinen Voraussetzungen finden kann. 497 


KLEINE MITTEILUNGEN 


Das Aufheizen einer Wand durch kon- Da die linke Seite dieser Identität von der Zeit 
stante Wärmestromdichten. 1. Die Lö- unabhängig ist, muß 
sungsmethode. Die Lösung dieser Aufgabe, 
die Zd. Lastovkat) zeichnerisch behandelt hat, ak, 
gelingt rechnerisch nach dem Summenansatz da- dr? 
durch, daß die Beharrungstemperatur als lineare 
Funktion der Zeit angenommen wird, deren Koeffi- 
zienten Polynome 3. bzw. 1. Grades der Raum- a? % 
koordinate sind. Der konstante Symmand in der jgentität C-° eine lineare Funktion, also k, ein 
Beharrungstemperatur und die Amplituden in der 2°’ 
Ausgleichstemperatur werden als Fourierkoeff-r Polynom 3. Grades der Raumkoordinate. Damit 
zienten bestimmt. wird aus der Annahme (1) 


2. Formulierung der Aufgabe. Für die Op—khotKroctkgnr?+t Kor? t (kat kt 
Wärmeleitung gilt die Differentialgleichung 


alsu der Koeffizient k, eine lineare Funktion 
der Raumkoordinate sein. Dann ist nach der 


u Dieses Polynom wird in die Differentialgleichung 
0» u 9 eingesetzt: 
dt 2° 
kathır=aldkunt6kur). 
Die Wanddicke wird mit ö bezeichnet. Rand- " E h 
bedingungen: Durch Koeffizientenvergleich folgt 


Ki = Eok.. Ku —ba ku. 


(52) —, 5 (22) ==q, 
Or/x=U Or/x=d Diese Werte werden in das Polynom eingesetzt: 
Anfangsbedingung: (d),., =. ISp—hotKkort Kzt?t Kyd? \ 


3. Berechnung eines partikulären Inte- + 2a (ka t3ku tl" 
grals. Da die Wand durch beide Begrenzungs- 
ebenen ständig aufgeheizt und keine Wärme ab- Dieses Polynom wird in die Randbedingungen ein- 
gegeben wird, muß auch die Beharrungstempe- gesetzt: 
ratur im Laufe der Zeit wachsen. Die einfachste ® q 
ot 6ukut=—- 


Annahme ‚ist die, daß die Beharrungstemperatur 
eine lineare Funktion der Zeit ist, deren Koeffi- 
zienten zunächst beliebige Funktionen der Raum- kot2dknt 35° hy + 6a Ayut 
koordinate sind: 


H=ekhd+hDdt...... WM. 


Dieser Ansatz wird in die Differentialgleichung ko= 
eingesetzt: 


Fa 
Durch Koeffizientenvergleich folgt 


dtr, , dek 
kı= al‘ te ri} !). Diese Werte werden in Gl. (2) eingesetzt: 

1), Zd. LaStovka: Forsch. Ing.-Wes. Bd. 13 (1942), I ku + 4 ( - + £ + 9 e m 
S, 262 bis 263. p 4 h) Ö 
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4. Berechnung des allgemeinen Integrals 
der Differentialgleichung und der Rand- 
bedingungen ohne Störungssummanden. 
Der Bernoullische Produkt- und d’Alembert- 
sche Exponentialansatz 

d,—e-“k’t(4sinkr+Bcoskz) 


[7] 


Störung 


s- 


wird in die Randbedingungen ohne 
summanden eingesetzt: 


A—=B, sinkö=0. 


Daher ist die Ausgleichstemperatur 


= n 2 
5, P I t » 
= Y B„e \F) cos aa 4 


— 
nl 


5. Erfüllung der Anfangsbedingüng. 


Die Temperatur 


Q 
dB—=kaot ; 


e | 
y —(ä - 
+ >) B, € 


a 


n==1 
wird in die Anfangsbedingung eingesetzt: 


ne 
> nnı _ q ( =*\ 


Koo + d } . 


1. B,. c08 3 Fl 
Mit der Orthogonalität folgt 

. d 
-5/6-5 


v 


9 
7 


6. Ergebnis. 


z- 11 
+ 


46 v 
RT ul n? 
n=2,4,6, +... 


-2+35 +23 


„ey 


u. or ral 


e 


Die Beharrungstemperatur 
G 
nm tlö 44284 


wird durch die Ordinaten einer Parabel gegeben, 
deren Minimum in der Mitte liegt, und die mit kon- 
stanter Geschwindigkeit wachsen. Die Oberflächen- 
temperatur ist 


' a 
urn fl 2 - 


Im Fall- einer Ziegelwand von der 
kcal 
mhGrad’ 


7. Beispiel. 
Dicke 0,38 m, der Wärmeleitfähigkeit 0,75 
ii ui kcal 
der spezifischen Wärme 0,22 > 
kg Grad 
1 705 ke; t 
[AU) 18 
m? 


und dem spezi- 
fischen Gewicht die konstante Ge- 


schwindigkeit, mit der die Ordinaten anwachsen, 
2q m? Grad 

1705-038 keal 

4 m?Grad 

1 

ER keal 


02. 
— 0,014 0 


und die größte Zeitkonstante 


9 375,1: 0,38° 


u h=1h49m4, s 


&. Lösung nach der Operatorenmethode. 
Die partielle Differentialgleic hung mit der Anfangs- 
bedingung geht durch die La plac esche Integral- 
transformation 


> 
d—=le-':tddt 


vu 


Re(s) >0 


in die gewöhnliche Differentialgleichung 
ER ;.. 


sv=a 
dx? 


über. Deren allgemeines Integral ist 


9 —K, Sin | - 2+K,Coi | 


Die transformierten Randbedingungen 


a 
i EE oe 8 


. (d Ö\ 
12), = - ar 


werden durch die transformierte Temperatur 


Si &(e 


a 


S au: s 
Sin I/ . 
az 


erfüllt. Die Temperatur wird nach der Riemann- 
Mellinschen Umkehrformel 


o+wi 
P 1 f ” 
= 5—, | etsJds 
Te 
o— @i 


berechnet; hierin muß der Integrationsweg so ge- 
wählt werden, daß die Pole der transformierten 
Temperatur links vom Integrationswege liegen. 


FOR s 
Sin 
a 


Die Pole der transformierten Temperatur liegen 
' „3:8 
in den Punkten s=0 unds=—? je Daher 
wird als Integrationsweg die Gerade s—=o > ge- 
wählt. Dieser Integrationsweg wird zunächst 
durch 2 Punkte symmetrisch zur reellen Achse be- 
grenzt. Diese beiden Punkte werden durch einen 
Bogen des Kreises miteinander verbunden, dessen 
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Gewählt wird der 


"Mittelpunkt der 0-Punkt ist. 
Auf 


Bogen, für dessen Punkte #e(s)<o ist. 
diesem Kreisbogen ist 


ets Sn et, 


Die Polarkoordinaten der Integrationsvariablen 
werden mit r und @ bezeichnet. 


co); [? _ 5) 


r r( 6\ 
2 008 5 ii , 


r gpÖ 
nV Es? 


gekürzt, so ist der absolute Betrag des Zählers S1 
und der des Nenners 


RE 7 
cos? V: sin . ; 
MRERGEEN Erz 


Sof V - _ c08 % 


/ 
/ 


Läßt man den Kreisbogenhalbmesser die Folge 


(2n — 1)? n? Fr durchlaufen, so ist 


Pi: op 6 n.o 1 
si/ _, . 2 (9n — 22 @ Bi... 
cos | „sin, > cos? (2n )yzsinzsz; 
in einem gewissen Intervall mit dem Mittel- 
punkte g=n. Die untere Grenze 9, dieses Inter- 


valls ist durch die im 1. Quadranten liegende 
Wurzel En der Gleichung 


gegeben. In jenem Intervall , Sp sS2n — gq, ist 


nn 


/ cos? V 5 sin » 
N A ee 2V;- 


D 


Z— 


Außerhalb jenes Intervalles ist 


’ er 
cos? | sin — 
a 2 


Cof? | 


> Tg l ä co 
a 


Tg (2n 1) 2 l 1— sin? In 


nn] 


| | — Tg V 


- Tg | 7 er — )n— z ya 


4 


Daher geht das Integral über den Kreisbogen 
nach 0. Mithin kann das Integral über die Strecke 
durch das Integral über den Kreisbogen zu einem 
Integral über eine geschlossene Kurve ergänzt 
werden. Da der Integrand eine eindeutige Funktion 
der Integrationsvariablen ist, kann das Integral 
über die geschlossene Kurve nach dem Residuen- 
satz ausgewertet werden. Um das Residuum für den 
Pol s=0 zu bestimmen, werden die Exponential- 
funktion und die Hyperbelfunktionen nach Potenzen 
der Integrationsvariablen entwickelt. 


Das Residuum für den Pols=0i 


a x? ö 
>) nn 
25:75 +4: 


Für die Pole s= — j? n? m sind die Resıduen 


n?a 


6 ö 
amatcos iz (z- 5) 
Pin Wirken 2 2143(- 1 


‚nd Ö 
cj55 


jax 
5; 


2 at 
4 -Pez 


>» ? Pr öde cos 


/ > 
r - 
Coj? V — (8 — — 
- s 
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